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ABSTRACT 

This thesis takes five FORTRAN IV programs from "Com- 
puter Programs for Computational Assistance in the Study 
of Linear Control Theory" by Melsa and Jones and translates 
them into a microcomputer BASIC language to run on an inex- 
pensive microcomputer system. Three of the five programs are 
state variable programs. They are BASMAT for basic matrix 
manipulation, RTRESP for rational time response, and GTRESP 
for graphical time response. Two are transfer function pro- 
grams, FRESP for frequency response and RTLOC for root locus. 
A user's guide and example are included for each. 

A final example is used to demonstrate the utility 
the two transfer function programs as an aid to direct 
digital design in the w' -plane. 
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I. 



INTRODUCTION 



The purpose of this thesis is two fold. First, it 
is to show that standard computer programs useful in the study 
of linear control theory may be adapted to run on an inex- 
pensive home /microcomputer . To demonstrate this, five 
programs were chosen from [Ref. 1], three state variable 
programs and two transfer function programs. The three state 
variable programs are BASMAT, a basic matrix manipulation 
program, RTRESP, a rational time response program, and GTRESP, 
a graphical time response program. The two transfer function 
programs are FRESP, a frequency response program, and RTLOC, 
a root locus program. These programs as they appear in [Ref. 
1] are written in basic FORTRAN IV language to be run on 
a main frame computer system utilizing standard graphics 
subroutines. These programs were modified and rewritten in 
a microcomputer BASIC language which is an interpreted 
language. Generally these programs are limited to systems 
of 10th order. It is felt that this limitation is more than 
acceptable for the purposes of this thesis. Also no major 
effort has been made to analyze or improve the efficiency of 
the numerical methods used. An effort of this type is 
advised if these programs are to be modified for higher 
order systems. 
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The second part of this thesis investigates the methods 
and relationships involved in direct digital design in the 
w' -plane. The two transfer function programs adapted to run 
on a microcomputer from the first part of this thesis, FRESP 
and RTLOC , are used to aid in this investigation. 

Section II summarizes the common problems and consi- 
derations involved in the translation of the five programs 
to be run on a microcomputer system. Section III describes 
the three state variable programs, BASMAT, RTRESP, and 
GTRESP, and gives an example of their use and output. 

Section IV is similar to section III and describes the two 
transfer function programs, FRESP and RTLOC. 

Section V deals with the w'-plane. Subsection A gives 
some background on the w'-plane and subsection B provides a 
simple example using the two transfer function programs to 
compare the s and w' planes. Subsection C develops templates 
of some constant parameters in the w'-plane. Subsection D 
ends the section with a more involved example. 

Section VI provides some conclusions and recommendations. 
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II. TRANSLATION CONSIDERATIONS 



The five programs in this thesis, BASMAT, RTRESP, GTRESP, 
FRESP, and RTLOC , were translated directly from programs of 
the same name in [Ref. 1]. These programs were written in 
the basic FORTRAN IV computer language, a compiled language 
designed to run on a main frame or mini computer system. In 
the first part of this thesis these programs are translated 
into a microcomputer BASIC language adapted to run on the 
microcomputer system described in Appendix A. 

In general the translation from FORTRAN IV to the micro- 
computer BASIC used posed no serious problems. Even though 
the BASIC language used is by necessity an abridged version 
of the BASIC language found on larger more expensive systems, 
it was extensive enough to provide the necessary commands 
for these programs. 

The input portion of the programs were rewritten to be 
interactive for convenient keyboard entry of problem para- 
meters eliminating input formatting errors. The programs 
requiring more extensive input were given the added capability 
of saving and retrieving a problem description from a disk 
file for the user's convenience. 

Unlike FORTRAN IV the BASIC used does not have complex 
math capabilities. Therefore, mathematical operations on 
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complex quantities were programed separately for the real and 
imaginary parts. 

Another important difference in the languages affecting 
translation is that BASIC has no provision for local variables 
All variables in the main program and all subroutines are 
global. This caused some bookkeeping problems in translating 
subroutines to prevent undesired side effects. Another compli 
cation was that FORTRAN IV considers the first four characters 
of a variable name for identification of that unique variable. 
The BASIC used considers only the first two. This further 
complicated the bookkeeping of variables and resulted in 
variable names being assigned just because they were different 
from the rest and with no relation to the quantity represented 

The output routines were written to conform as close as 
possible to the FORTRAN IV version. Since the microcomputer 
system used is limited to eighty columns on eight and a half 
inch wide paper, provisions were made within the program to 
automatically switch to a condensed character font when 
necessary to output greater than eighty characters per line. 

The programs requiring graphical output created some 
unique problems. First, no standard library subroutines for 
graphics as used in the FORTRAN IV programs were available 
on the microcomputer system used. As a result all graphics 
routines had to be developed for the system used. Also the 
appropriate variables were redimensioned giving consideration 
to the resolution of the graphics available on the 
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microcomputer system to optimize the program somewhat. 
Secondly, although it is possible to mix graphics and text 
on the microcomputer system, it is not done in a straight for- 
ward manner and utilizes extra memory. For memory considera- 
tions and convenience it was decided to provide sufficient 
information for interpretation of the graphical output 
below each graph. It is felt that this method is satisfac- 
tory and creates little inconvenience to the user. To 
enhance the interpretation of the graphical output much 
time and effort was devoted to developing plotting routines 
to display the data in relation to axes labeled with tic-marks 
and boundaries of integer vice fractional values. 

A major obstacle that was overcome was the identification 
and correction of a memory management problem unique to the 
microcomputer system used. This problem affected only those 
programs requiring graphical output. The essence of the 
problem was due to the size of the programs and the number 
of variables used, parts of the program and stored variables 
were being over written in memory by the graphics routines. 
This problem was finally solved by making appropriate changes 
in the memory management scheme. For a more detailed explana- 
tion of the solution see Appendix A. 
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Ill . 



STATE VARIABLE PROGRAMS 



The three state variable programs discussed in this sec- 
tion are modified versions of the programs of the same names 
found in [Ref. 1]. These programs may be used as tools for 
the analysis and design of linear control systems expressed 
in the following state variable form: 

x ( t ) = Ax ( t ) + bu ( t ) 
u(t) - K[r(t) - k^x ( t ) ] 
y ( t ) = c T x(t) 

The Basic Matrix program (BASMAT) described in part A 
of this section is used to compute the determinant, inverse, 
characteristic polynomial, and eigenvalues of the square 
matrix A. It is also used to determine the resolvent matrix 
(sI-A) ^ and state transition matrix exp(At). 

The Rational Time Response program (RTRESP), described 
in part B, is used to determine a closed form expression 
for the time response of a system. The input function r(t) 
is required to have a rational time response and no repeated 
eigenvalues are allowed in the combination of the system 
and input . 

The Graphical Time Response program (RTRESP), described 
in part C, is used to produce a graphical display of the time 
response of a system to an arbitrary input. 
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RTRESP and GTRESP can be used to study open loop systems 
by letting K equal zero and unforced systems by letting r(t) 
equal zero. 
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A. BASMAT 



1 . Basic Matrix Program (BASMAT) 

a. Introduction 

When done by hand, matrix manipulations can be 
quite tedious and the chances of an error being made are 
great. In the study of linear state variable analysis, a 
computer program to do these manipulations is almost 
essential. It can do the necessary manipulations much more 
quickly and accurately and allow the user to devote his 
time and energy to design and analysis. 

b. Description of Program 

BASMAT [Ref. 1: pp . 7,8] takes a matrix A, 

and computes the determinant of A (det A), the inverse of 
A (A ^ ) , the characteristic polynomial (det ( sj[-A) ) , and 
eigenvalues of A (X { ‘)as well as the resolvent matrix 
= (sI-A) -1 

and the state transition matrix 
<j6(t) = exp (At) 

The state transition matrix is expressed as matrix coeffi- 
cients times the natural modes exp(\t) with the complex 
eigenvalues expressed as damped sine and cosine terms. The 
resolvent matrix is written as 
<j6(s) = adj (sI-A) / det (si- A) 

and the matrix numerator, adj(s < I-A), is output as matrix 

coefficients of powers of s so that it takes the form 

adj(sI-A) = F- + F 0 s + ... + F„ T s^ - ^ 

~ — 1 Z N 



14 



The BASMAT program interactively accepts input 



of a matrix and calls subroutines to perform the appro- 
priate calculations. The subroutines used are CHREQ, CHREQA, 
PROOT , DET, SIMEQ, and STMST . These subroutines are listed 
below with a brief description. 

CHREQ. This subroutine is used to determine the 
characteristic polynomial det(sI-A), and the resolvent 
matrix (s.I-A) - ^ for the matrix A. The Leverrier algorithm 
is used to compute the resolvent matrix 



where D(s) is the characteristic polynomial det(sI-A). 

The coefficients of the characteristic polynomial 
are determined by subroutine CHREQA. 

CHREQA. This subroutine is called by CHREQ to 
determine the characteristic polynomial, det(sI-A), for the 
A matrix. The principal-minor method is used. 

PROOT. This subroutine uses a modified Bair- 
stow method for root extraction to determine the roots of a 
polynomial with real coefficients. It is used to calculate 
the eigenvalues of the A matrix. 

DET. This subroutine computes the determinant 
of a matrix. A gauss elimination method to place the matrix 




in the form 



N 
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in upper triangular form is used. It is called by CHREQA to 
calculate sub-determinants. 

SIMEQ. This subroutine is used in the FORTRAN 
version of BASMAT in [Ref. 1: pp. 7,8] to determine the 

inverse of the A matrix. In the BASIC version a machine 
language subroutine [Ref. 2] is substituted for speed 
and convenience. 

STMST. This subroutine is used to compute the 
state transition matrix 
<^>(t) = exp (At) 

for a matrix A. It uses the Sylvester Expansion Theorem. 
Eigenvalues of the A matrix must be provided. This routine 
can not handle duplicate eigenvalues. 

c. Program Translation Problems 

See section II, Translation Considerations. 

2 . BASMAT User's Guide 

This program occupies about 7k bytes of memory and 
does not utilize the graphics pages, therefore it is not 
necessary to relocate the disk operating system. 

Push "return" after each input. 

STEP 1: 

SCREEN PROMPT: 

BASIC MATRIX PROGRAM 
PROBLEM IDENTIFICATION -> 
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REMARKS : 



Label the problem. Type in any appropriate combina- 
tion of letters, numbers, and/or symbols, excluding 
commas. This input is limited to 255 characters. 

STEP 2: 

SCREEN PROMPT: 

PLANT ORDER? -> 

REMARKS : 

Enter the number that corresponds to the order of 
the A matrix. Maximum order allowed is ten. 

STEP 3: 

SCREEN PROMPT: 

INPUT PLANT MATRIX. 

A(row,column)= 

REMARKS : 

Input the A matrix by typing in each element as 
prompted. The program will ask for each element 
beginning with the first row and going from left to 
right . 

STEP 4: 

SCREEN PROMPT: 

HARDCOPY? (Y/N) -> 

REMARKS : 

After this step is completed the program will output 
the A matrix for reference. This will be followed by 
the inverse of the A matrix, the determinant of the A 
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matrix, the matrix coefficients of the resolvent 
matrix numerator, the characteristic polynomial 
coefficients, the eigenvalues of the plant matrix, 
and finally the elements of the state transition 
matrix. See Figure 1 for a sample output. 

3 . BASMAT Example 

This example will use the system matrix 



A 




1 0 
0 1 
-3 -3 



to demonstrate the use of the BASMAT program. It will refer 
to the steps described in the previous section, BASMAT 
User's Guide. 

Step 1. Enter "EXAMPLE <CR>" where <CR> = return. 

Step 2. Enter "3 <CR>" 

Step 3. Enter "0 <CR> 1 <CR> 0 <CR> 0 <CR> 0 <CR> 1 <CR> -2 
<CR> -3 <CR> -3 <CR> " 

Step 4. Enter "Y <CR>" 

The resulting output is shown in Figure 1. These 
results are interpreted as follows: 

The plant matrix is shown just as it was entered with 
rows horizontal and columns vertical. This is the format 
used for all matrix output. 

The inverse of the plant matrix is shown followed by 
the scalar determinant of the A matrix. 
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The matrix coefficients of the resolvent matrix 



numerator are given as powers of s. The characteristic 
polynomial is listed as coefficients of powers of s. 

From this output , the resolvent matrix 
<^>(s) = adj (sI-A) /det (sI-A) 

may be written. The characteristic polynomial, det(sI-A), 

2 3 

is 2+3s+3s +s . The first element of the resolvent matrix 

2 

numerator is 3+3s+s making the first element of the 
resolvent matrix 

2 

As / x 3+3s+s 

= 2 3 

; 2+3s+3s Z +s J 

The real and imaginary parts of the eigenvalues of 
the plant matrix are listed. Finally, the elements of the 
state transition matrix are given. The first element of 
the state transition matrix can be written 

<^>(t) = 0.333e“ 2t + 0.667e"°* 5t cos 0 . 866t+l . 15e“° ,5t sin 0.866t 
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BASIC MATRIX PROGRAM 

PROBLEM I DENT I F I CAT I ON-EXAMPLE 



XXXXXX X X X X X X X X X X X X X X X X X X X X X X X X X X X 

THE PLANT MATRIX 



0 

0 

_o 



1 

0 



0 

1 



X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X 

THE INVERSE OF THE PLANT MATRIX 



-1.5 

1 

-O 



-1.5 

0 

1 



-.5 

0 

0 



X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X 



THE DETERMINANT OF THE PLANT MATRIX 



X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X 

THE MATRIX COEFFICIENTS OF THE RESOLVENT MATRIX 



THE MATRIX COEFFICIENT OF 3" 2 



1 

0 

0 



0 

1 

0 



0 

0 

1 



THE MATRIX COEFFICIENT OF S A 1 
1 

-3 

THE MATRIX COEFFICIENT OF S A 0 



3 

0 

-2 



0 

1 

0 



o 

-2 

0 



3 

0 

_o 



1 

0 

0 



X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X 
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NUMERATOR 



Figure 1 



BASMAT Output 



Figure 1 (Cont.) 

THE CHARACTERISTIC POLYNOMIAL IN ASCENDING POWERS 



2 3 3 1 

X X * X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X 



THE EIGENVALUES OF THE PLANT MATRIX 
REAL IMAGINARY 

-2 0 

-.5 -.866025464 

-.5 .366025404 



X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X 



X 



THE ELEMENTS OF THE STATE TRANSITION MATRIX 



THE MATRIX COEFFICIENT OF EXP<-2)T 
.333333333 .333333333 

— . o o 6 6* 6' 6 6 6 — • $ c> 6 o o o & 6 

1.33333333 1.33333333 



~ ■ Ol , v ' ju 666 u 

1.33333333 



THE MATRIX COEFFICIENT OF EXP< - . 5) TXCOS< . 366025404) T 
. 666666666 - . 333333333 - . 333333333 



• 6<u?6£66oo6 



i < Ouu’JuCu / 



6 uOOuC'ju 



THE MATRIX COEFFICIENT OF EXP < - . 5) TX 3 1 N < . 8666 2546 4> T 



1.15476054 1.73285081 

-1.15476054 -.577356268 

9 -1.15470854 



57735 



8263 



- . 577350268 



¥ V X * X X X X X XX X X X X X X V X X X X X X X X X X XX XXX 

. ^ . \ A A A /\ .rs a *\ .■% .a ✓ % ,s A A /s A A /•. A ,*% A A ✓% A 



Figure 1 
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B. RTRESP 



1 . Rational Time Response Program (RTRESP) 

a. Introduction 

Frequently it is desirable to know the response 
of a system as a function of time. A computer program can 
determine this quicker and more accurately than by hand. 

b. Description of Program 

RTRESP determines the time response in closed 
form of the closed loop system 

x ( t ) =Ax (t)+bu(t) 
u(t)=K[r(t) -k T x ( t ) ] 
y (t )=c T x(t ) 

due to any initial conditions x(0) and input r(t) for t>0. 

The system must have a rational Laplace transform R(s) with 

a pole-zero excess of at least one. [Ref. 1: pp. 11,12] 

"t h 

The input r(t) is treated by forming a m -order 
dynamic system whose initial condition response is equal to 
r(t) for a specific set of initial conditions. This system 
is combined with the original system and then the complete 
response in closed form is determined from the subroutine 
STMST. The order of the combined system must be ten or 
less . 

Various primary and utility subroutines are used 
in RTRESP. The primary subroutines are listed below with a 
brief description: 
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CHREQA. This subroutine determines the charac- 
teristic polynomial, det(sI-A), for the matrix A using the 
principal-minor method. 

PROOT. This subroutine uses a modified Bairstow 
method for root extraction to determine the roots of a poly- 
nomial with real coefficients. It is used to determine the 
roots of the numerator and denominator of R(s) when they are 
entered in polynomial form. PROOT is also used to determine 
the eigenvalues of the combined system matrix (i.e. the roots 
of the characteristic polynomial determined by CHREQA). 

SEMBL. This subroutine determines the coeffi- 
cients of a polynomial from the roots of the polynomial. 

It is used to determine the polynomial coefficients of the 
numerator and denominator of R(s) when they are entered 
in the factored form. This subroutine together with PROOT 
provide the feature that R(s) may be entered in either 
factored or polynomial form. It will appear in the output 
in both forms. 

STMST. This subroutine computes the state transi- 
tion matrix 
(£>( t) = exp(At) 

for a matrix A. It uses the Sylvester Expansion Theorem. 
Eigenvalues of the A matrix must be provided. This routine 
can not handle duplicate eigenvalues, therefore it is 
necessary that the combined system and input have no 
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repeated eigenvalues. STMST is used to determine the state 
transition matrix of the augmented system. 

c. Program Translation Problems 

See section II, Translation Considerations. 

2 . RTRESP User's Guide 

Push "return" after each input. 

STEP 1: 

SCREEN PROMPT: 

RATIONAL TIME RESPONSE PROGRAM (RTRESP) 

THIS PROGRAM DETERMINES THE TIME RESPONSE OF A 
CLOSED-LOOP SYSTEM DUE TO SPECIFIED INITIAL 
CONDITIONS AND INPUT. SYSTEM MUST HAVE A RA- 
TIONAL LAPLACE TRANSFORM WITH A POLE-ZERO EXCESS 
OF AT LEAST ONE. 

PROBLEM IDENTIFICATION -> 

REMARKS : 

Label the problem. Type in any appropriate combina- 
tion of letters, numbers, and/or symbols, excluding 
commas and colons. This input is limited to 255 
characters. This is also the name used to save and 
retrieve the disk file containing the problem 
description. (See steps 2 and 3) 

STEP 2: 

SCREEN PROMPT: 

WAS THIS PROBLEM DESCRIPTION PREVIOUSLY 
SAVED? (Y/N) -> 

REMARKS : 

This step provides the option of recalling a pre- 
viously saved problem description from the disk. 
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If "Y" is typed, the program looks for the problem 
description saved under the label typed in STEP 1. 
The program then runs to completion using the 
retrieved program description. If an "N" is 
typed, the program continues with STEP 3. 

STEP 3: 

SCREEN PROMPT: 

DO YOU WANT TO SAVE THIS PROBLEM DESCRIPTION? 
(Y/N) -> 

REMARKS : 

A positive response will save the problem descrip- 
tion, that will be entered in the following steps, 
to the disk under the label entered in STEP 1. 

This problem may be retrieved during a later session 
by entering the proper label in STEP 1 and typing 
"Y" in STEP 2. 

STEP 4: 

SCREEN PROMPT: 

ORDER OF THE SYSTEM? -> 

REMARKS : 

The total order of the system and the input must 
not exceed ten. 

STEP 5: 

SCREEN PROMPT: 

INPUT SYSTEM (A) MATRIX 
A( 1 , 1 ) = 
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REMARKS : 



Enter the elements of the A matrix as prompted. 

The format is A(row , column) = . 

STEP 6: 

SCREEN PROMPT: 

THE A MATRIX 

(display of the A matrix) 

ANY CHANGES? (Y/N) -> 

REMARKS : 

If an "N" is typed, the program proceeds to STEP 9. 
To correct the matrix type "Y" and the program will 
proceed to STEP 7. 

STEP 7: 

SCREEN PROMPT: 

TYPE ROW, COLUMN OF THE ELEMENT TO BE 
CORRECTED -> 

STEP 8: 

SCREEN PROMPT: 

A (row, column) = 

REMARKS : 

Enter the correct value. The program will return 
to STEP 6. 

STEP 9: 

SCREEN PROMPT: 

INPUT THE CONTROL (B) VECTOR 
B ( 1 ) = 
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REMARKS : 



After the B vector is entered, a correction sequence 
similar to steps 6 through 8 is encountered. 

STEP 10: 

SCREEN PROMPT: 

INPUT THE OUTPUT (C) VECTOR 
C(l) = 

REMARKS : 

After the C vector is entered, a correction sequence 
similar to steps 6 through 8 is encountered. 

STEP 11: 

SCREEN PROMPT: 

INPUT THE FEEDBACK COEFFICIENTS 
FEEDBACK COEFFICIENT (1) = 

REMARKS : 

After the feedback coefficients are entered, a 
correction sequence similar to steps 6 through 8 
is encountered. 

STEP 12: 

SCREEN PROMPT: 

INPUT THE GAIN -> 

REMARKS : 

Enter the controller gain. 
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STEP 13 



SCREEN 


PROMPT : 




INPUT THE INITIAL CONDITIONS 
X0(1) = 



REMARKS : 

After the initial conditions are entered, a 
correction sequence similar to steps 6 through 
8 is encountered. 



STEP 14: 




SCREEN 


PROMPT : 

ENTER THE INPUT GAIN -> 



REMARKS : 

Enter the gain of the input function, R(s). 



STEP 15: 




SCREEN 


PROMPT : 

'KEY' = P FOR POLYNOMIAL COEFFICIENT FORM 
'KEY' = F FOR POLYNOMIAL FACTORED FORM 

INPUT 'KEY', ORDER FOR NUMERATOR POLYNOMIAL -> 



REMARKS : 

Choose the preferred method of entering the numerator 
polynomial of the input function, R(s). 

STEP 16 A: 



SCREEN 


PROMPT (assumes coefficient form chosen): 
INPUT POLYNOMIAL COEFFICIENTS 




INPUT COEFFICIENTS IN ASCENDING PWRS OF S 
INPUT COEFFICIENT OF S^O -> 
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REMARKS : 



Enter the coefficient along with its algebraic 

sign as prompted. The program assumes that the 

coefficient of the highest power of s is one. 

STEP 16B (assumes factored form chosen): 

SCREEN PROMPT: 

INPUT POLYNOMIAL FACTORS 

FOR FACTOR 1 
REAL PART -> 

IMAGINARY PART -> 

REMARKS : 

Enter the factor with its algebraic sign. Do not 
enter the associated root. (e.g. For "(s-1)" enter 
"-1".) After the real part of the factor is entered, 
the program asks for the imaginary part , allowing 
for input of quadradic factors. If a non-zero 
imaginary part is entered, the program automatically 
enters its complex conjugate. 

STEP 17: 

SCREEN PROMPT: 

•KEY' = P FOR POLYNOMIAL COEFFICIENT FORM 
'KEY' = F FOR POLYNOMIAL FACTORED FORM 

INPUT 'KEY', ORDER FOR DENOMINATOR POLYNOMIAL -> 

REMARKS : 

The denominator is entered just as described for 
the numerator in steps 15 and 16. 
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After this step is completed, the program will 
output all input for reference followed by the 
time response of x(t) and the system output, 
y( t ) . 

3. RTRESP Example 

It is desired to know the time response of the 
closed loop system 



x(t ) = 


" 0 


1 


x(t ) + 


"o" 


**. 


_-l 


-1 




1 



u ( t ) = 5{r( t ) - [1.0 0.5] x ( t ) } 

y(t) = [1 1] x ( t ) 

if the input function is given by 

R(s) = 2.0 (1/s) 
and the initial conditions are 

x(t ) =0 

Referring to the steps described in the pervious sec- 
tion, RTRESP User's Guide, this problem would be entered as 
follows : 

Step 1. Enter "EXAMPLE <CR>" where <CR> = carriage 

return . 

Step 2 and Step 3. Enter "N <CR> . " 

Step 4. Enter "2 <CR>." 

Step 5. Enter "0 <CR> 1 <CR> -1 <CR> -1 <CR>." 

Step 6. Enter "N <CR> . " 
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Step 


9. 


Enter " 


'0 <CR> 1 <CR> . " 


Enter 


appropriate 


response for 


correction 


sequence . 






Step 


10. 


Enter 


"1 <CR> 1 <CR> . " 


Enter 


appropriate 


response for 


correction 


sequence . 






Step 


11. 


Enter 


"1 <CR> .5 <CR> . 


" Enter appropriate 


response for 


correction 


sequence . 






Step 


12. 


Enter 


"5 <CR> . " 






Step 


13. 


Enter 


"0 <CR> 0 <CR> . " 


Enter 


appropriate 


response for 


correction 


sequence . 






Step 


14. 


Enter 


"2 <CR> . " 






Step 


15. 


Enter 


"P,0 <CR> . " 






Step 


16A . 


Enter 


"1 <CR> . " 






Step 


17. 


Enter 


"F , 1 <CR> 0 <CR> 


M 

• 





The program will now run to completion resulting in the out- 
put shown in Figure 2. 

The output begins with the program name and problem 
identification. The A, B, and £ matrices are shown for 
reference followed by the feedback coefficients and the 
controller gain. Next, the initial conditions and input 
gain are listed. The numerator and denominator polynomials 
of the input function are given in both polynomial coeffi- 
cient form and factored form regardless of the method in 
which they were entered. The polynomial coefficient form is 
given as a list of coefficients from left to right in 
ascending powers of s. The highest power of s is always 
one. The polynomial factored form appears as a list of the 



31 



real and imaginary parts of each root of the polynomial. 

Notice that the numerator is only given in coefficient 
form since it is of zero order and no roots exist . 

The time response of the state x(t) is given in the 
form of vector coefficients of the various natural modes of 
the system and the input function. From this output 
x(t) is seen to be 

— 1 TSt —1 7 R-f- 

X x (t) = -1 . 67e ,/OL cos 1 . 71t - 1.70e L s i n 1.71t + 1.67 

-1 7 c it 

x 2 (t) = 5.83e * sin 1.71t 

The time response of the output y(t) is given as scalor 
coefficients of the same natural modes as x(t) so it is seen 
that 

y ( t ) = -1.67e" 1,75t cos 1.71t + 4.13e _1 * 75t sin 1 . 71t + 1.67 
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RATIONAL TIME RESPONSE 

PROBLEM IDENTIFICATION -> EXAMPLE 

************ * * * * ****** * * * * * * * * * * * * * * * * * * * * * * * * * * * * 

THE A MATRIX 

8 1 

-1 -1 

THE B MATRIX 

0 1 

THE C MATRIX 

1 1 

FEEDBACK COEFFICIENTS 
1 .5 

GAIN = 5 

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 
INITIAL CONDITIONS, XC0) 

0 0 

RGAIN = 2 



NUMERATOR POLYNOMIAL OF PCS) - ASCENDING F'WRS OF S 
1 



Figure 2 RTRESP Output 
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Figure 2 (Cont.) 



DENOMINATOR POLYNOMIAL OF R<S) - ASCENDING PWRS OF 
9 1 

DENOMINATOR ROOTS ARE 

REAL PART IMAGINARY PART 

0 0 

X X X X X X X X X X X X X X X X X X XX XX X X X X X X X X X X X X X X X X XXX X X X X X X X X X 

THE TIME RESPONSE OF THE STATE X<T) 



THE VECTOR COEFFICIENT OF EXP< - 1 . 75) TSCOS< 1 . 7 139 13c 

— 1 • 66666667 0 

THE VECTOR COEFF I C I ENT OF EXP < - 1 . 75) TX S I N < 1 . 7 1 39 1 3c 
-1.70175824 5.83459966 

THE VECTOR COEFFICIENT OF EXP<0>T 
1.66666667 9 . 3 1 322575E- 1 0 

X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X 

THE TIME RESPONSE OF THE OUTPUT Y(T) 

THE COEFFICIENT OF EXP < - 1 . 75) TX COS < 1 .71391 365) T 

- 1 . 66666667 

THE COEFF I C I ENT OF EXP < - 1 . 75) T X S I N < 1 . 7 1 39 1 365) T 
4. 13284143 

THE COEFFICIENT OF EXPORT 
1 . 66666667 

* X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X 



Figure 2 



S3) T 
S5) T 
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C. GTRESP 



1 . Graphical Time Response Program (GTRESP) 

a. Introduction 

Knowing the response of a system as a function of 
time can be very helpful in the study and analysis of that 
system. The graphical display of this response can give 
much insight into a system's response characteristics. 

b. Description of Program 

GTRESP [Ref. 1: pp. 22-28] determines and 

graphically displays the time response of the closed loop 
system 



x(t) = Ax(t) + bu(t) 
u(t) = K [r(t) - k T x ( t ) ] 
y(t) = c T x(t) 

due to any initial conditions x(0) and input r(t). 

The basic purpose of this program is similar 
to that of the RTRESP program described earlier. The 
difference between the two programs is that GTRESP will 
determine the time response for arbitrary input functions 
which may not have rational Laplace transforms but RTRESP 
requires a rational Laplace transform. Also GTRESP produces 
a graphical display of the time response instead of a 
closed form expression as in the RTRESP program. 

The GTRESP program uses a fourth-order Runge-Kutta 
numerical integration algorithm to calculate the time response. 
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The main subroutines used are CALCU , RUNGE, TRESP, and YDOT . 
Also various utility and plotting subroutines are used. 

CALCU. This subroutine is called by TRESP to 
determine the reference input r(t) and the control input 

u(t)=K [r(t) - k T x(t)] 

The reference input r(t) must be defined by the user by 
inserting the appropriate BASIC coding into the subroutine 
between line numbers 5010 and 5500. The reference input is 
represented by the variable R and the control input by the 
variable U. 

RUNGE. This subroutine is called by TRESP and 
contains the actual fourth-order Runge-Kutta integration 
algorithm. It must be executed four times for each inte- 
gration step. 

TRESP. This subroutine is the driving subroutine 
which calls the subroutines CALCU, RUNGE, and YDOT along with 
the necessary plotting routines. It calculates the time 
response of the closed loop linear system described by the 
input parameters and plots the desired variables. 

YDOT. This subroutine is called by TRESP to 
compute the derivative 
x(t )=Ax(t )+bu(t ) 

It is designed to handle linear systems but can easily be 
modified to handle nonlinear and time varying systems which 
would give GTRESP a nonlinear and time varying capability. 
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c. Program Translation Problems 

See section II, Translation Considerations. 

2 . GTRESP User's Guide 

Push "return" after each input. 

STEP 1: 

SCREEN PROMPT: 

GRAPHICAL TIME RESPONSE PROGRAM (GTRESP) 

THIS PROGRAM DETERMINES AND GRAPHICALLY 
DISPLAYS THE TIME RESPONSE OF A CLOSED- 
LOOP SYSTEM DUE TO SPECIFIED INITIAL 
CONDITIONS AND INPUT. 

PROBLEM IDENTIFICATION -> 

REMARKS : 

Label the problem. Type in any appropriate 
combination of letters, numbers, and/or symbols, 
excluding commas and colons. This input is 
limited to 255 characters. This is also the name 
used to save and retrieve the disk file containing 
the problem description. (See steps 2 and 3) 

STEP 2: 

SCREEN PROMPT: 

WAS THIS PROBLEM DESCRIPTION PREVIOUSLY 
SAVED? (Y/N) -> 

REMARKS : 

This step provides the option of recalling a pre- 
viously saved problem description from the disk. 

If "Y" is typed the program looks for the problem 
description saved under the label typed in STEP 1. 
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The program then runs to completion using the 
retrieved program description. If an "N" is 
typed the program continues with STEP 3. 

STEP 3: 

SCREEN PROMPT: 

DO YOU WANT TO SAVE THIS PROBLEM DESCRIPTION? 
(Y/N) -> 

REMARKS : 

The problem description will be entered in the 
following steps. 

A positive response will save the problem des- 
cription to the disk under the label entered in 
STEP 1. This problem may be retrieved during a 
later session by entering the proper label in STEP 
1 and typing "Y" in STEP 2. 

STEP 4: 

SCREEN PROMPT: 

ORDER OF THE SYSTEM? -> 

REMARKS : 

The total order of the system and the input must not 
exceed ten. 

STEP 5: 

SCREEN PROMPT: 

INPUT SYSTEM (A) MATRIX 
A(1 ,1) = 
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REMARKS : 



Enter the appropriate element of the A matrix. The 
format is A( row, column ) = . 

STEP 6: 

SCREEN PROMPT: 

THE A MATRIX 

(display of the A matrix) 

ANY CHANGES? (Y/N) -> 

REMARKS : 

If an "N" is typed, the program proceeds to STEP 9. 
To correct the matrix type "Y" and the program will 
proceed to STEP 7. 

STEP 7: 

SCREEN PROMPT: 

TYPE ROW, COLUMN OF THE ELEMENT TO BE 
CORRECTED -> 

STEP 8: 

SCREEN PROMPT: 

A (row, column ) = 

REMARKS : 

Enter the correct value. The program will return to 
STEP 6. 

STEP 9: 

SCREEN PROMPT: 

INPUT THE CONTROL (B) VECTOR 
B( 1 ) = 
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REMARKS 



After the B vector is entered a correction 
sequence similar to steps 6 through 8 is 
encountered . 

STEP 10: 

SCREEN PROMPT: 

INPUT THE OUTPUT (C) VECTOR 
C(l) = 

REMARKS : 



After the C vector is entered a correction sequence 
similar to steps 6 through 8 is encountered. 

STEP 11: 

SCREEN PROMPT: 

INPUT THE FEEDBACK COEFFICIENTS 
FEEDBACK COEFFICIENT (1) = 

REMARKS : 

After the feedback coefficients are entered a correc- 
tion sequence similar to steps 6 through 8 is 
encountered. 

STEP 12: 

SCREEN PROMPT: 

INPUT THE CONTROLLER GAIN -> 



REMARKS : 



Enter the controller gain 



STEP 13 



SCREEN PROMPT: 



INPUT 

X0(1) 

REMARKS : 


THE INITIAL CONDITIONS 


After the 


initial conditions are entered a correc- 



tion sequence similar to steps 6 through 8 is 
encountered. 

STEP 14: 

SCREEN PROMPT: 

ENTER THE FOLLOWING PARAMETERS: 

NOTE. . (TF-TZ) / (DT*FR) =< 100 

INITIAL TIME (TZ) “> 



REMARKS : 




Enter the 


initial time of the time interval of 


interest . 


Due to program constraints the initial 


and final 


times and the time step and frequency of 



output must be chosen so that they satisfy the 
relat ion . 

(TF-TZ) / (DT*FR) =< 100 

STEP 15: 

SCREEN PROMPT: 



FINAL 

REMARKS : 


TIME (TF) -> 


Enter the 
interest . 


final time of the time interval of 
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STEP 16: 



SCREEN PROMPT: 

TIME STEP (DT) -> 

REMARKS : 

Enter the time increment for each step. 
STEP 17: 

SCREEN PROMPT: 

FREQUENCY OF OUTPUT (FR) -> 



REMARKS : 

Enter an integer n and the program will print out 
data on every nth time step iteration. Ensure 
that the relation expressed in the remarks of STEP 
14 is satisfied. 

STEP 18: 

SCREEN PROMPT: 



YOU MAY 


PLOT UP TO 


8 VARIABLES 


VS TIME 


VARIABLE 


NUMBER 


VARIABLE 


NUMBER 


X1(T) 


1 


X8(T) 


8 


X2(T) 


2 


X9(T) 


9 


X3(T) 


3 


X10(T) 


10 


X4(T) 


4 


E(T) 


11 


X5(T) 


5 


U(T) 


12 


X6(T) 


6 


Y(T) 


13 


X7(T) 


7 


R(T) 


14 



HOW MANY VARIABLES TO PLOT? MAX=8 ~> 

REMARKS : 

You may plot up to 8 variables vs time. These 
variables will be plotted on the same plot in the 
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output. Enter the number of variables you want to 
appear on the plot. 

STEP 19: 

SCREEN PROMPT: 

TYPE THE VARIABLE NUMBER <CR> -> 

REMARKS : 

Type the number associated with the variable of 
interest. After carriage return <CR> is typed, the 
program will continue to prompt for another variable 
number until the number of variables that the user 
indicated in STEP 18 has been entered. 

STEP 20: 

SCREEN PROMPT: 

YOU CHOSE THE FOLLOWING VARIABLES: 

(list of variable numbers entered in STEP 19) 

DO YOU WANT TO MAKE ANY CHANGES? (Y/N) -> 

REMARKS : 

To make a change in the variables to be plotted 
type a "Y" and the program will return to STEP 
18. If an "N" is typed, the program will run to 
comp let ion . 

3 . GTRESP Example 

It is desired to determine and plot the time response 
of the error 

e(t )=r(t )-y (t ) 

the input r(t), and the state variable x 2 (t) for the system 
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~ 0 l" 


x(t) + 


'o' 


-1 -1 




1 



u(t) = 5{r(t) - [1.0 0.5] x(t)} 

y(t) = [1 1] x(t) 



The time interval of interest is 
0<t <4 

The iteration step size (DT) is chosen to be 0.01 and printed 
output is desired every ten steps so FREQ equals 10. The 
initial conditions are 
x(t )=0 

and the input function is 

r(t) = 5.0 if 0<t<l 

r(t) = 0 otherwise. 

In order to define this input function the following 

BASIC coding is inserted between line numbers 5010 and 5500 

in the subroutine CALC. 

5020 IF TZ > 1 GOTO 5040 
5030 R = 5: GOTO 5500 

5040 R = 0 

5500 REM END OF ROUTINE DESCRIBING R(T) 

Referring to the steps described in section 2, 

GTRESP User's Guide, this problem would be entered as 
follows : 

Step 1. Enter "EXAMPLE <CR>" where <CR> = carriage 

return . 

Step 2 and Step 3. Enter "N <CR>." 

Step 4. Enter "2 <CR>." 
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Step 


5. 


Enter 


"0 <CR> 1 <CR> -1 


<CR> - 


1 <CR> . " 


Step 


6 . 


Enter 


"N <CR> . " 






Step 


9. 


Enter 


"0 <CR> 1 <CR> . " 


Enter 


appropriate 


response for 


correction 


sequence . 






Step 


10. 


Enter 


"1 <CR> 1 <CR> . " 


Enter 


appropriate 


response for 


correction 


sequence . 






Step 


11. 


Enter 


"1 <CR> .5 <CR> . 


" Enter appropriate 


response for 


correction 


sequence . 






Step 


12. 


Enter 


"5 <CR> . " 






Step 


13. 


Enter 


"0 <CR> 0 <CR> . " 


Enter 


appropriate 


response for 


correction 


sequence . 






Step 


14. 


Enter 


"0 <CR> . " 






Step 


15. 


Enter 


"4 <CR> . " 






Step 


16 . 


Enter 


".01 <CR> . " 






Step 


17. 


Enter 


"10 <CR>." Ensure that 


the constraint 



given in step 14 of GTRESP User's Guide is met. That con- 
straint is (TF-TZ) / (DT*FR) = 100 

In this case 
( 4-0 ) / ( . 01*10 )=40 
so the constraint is met. 

Step 18. Enter "3 <CR>." 

Step 19. Enter "11 <CR> 14 <CR> 2 <CR> . " This causes 
the error e(t), the input r(t), and the state variable XgCt), 
respectively, to be plotted as desired from the problem 
statement . 
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The program will now run to completion resulting in 
the output shown in Figure 3. 

The output begins with the program name and problem 
identification. The A, B, and C matrices are shown for 
reference followed by the feedback coefficients and the 
controller gain. Next the initial conditions and time 
parameters are listed. 

The second page of output lists in tabular form the 
value of time t and the corresponding values of the output 
y(t), the control u(t), and all of the state variables. 

The user has no control over the variables output in this 
form. 

The third page of output is the graph itself. Below 
the graph is enough information to properly interpret the 
plotted data. 
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GRAPHICAL TIME RESPONSE 



PROBLEM IDENTIFICATION -> EXAMPLE 



* + * 
^ \ A / 



? ¥ ¥ : 



: * X 



£ ¥ ¥ ¥ ¥ ¥ ¥ : 



THE A MATRIX 



0 1 
-1 -1 



THE B MATRIX 
0 1 

THE C MATRIX 
1 1 

FEEDBACK COEFF I C I ENTS 
1 .5 

GAIN = 5 



INITIAL C OND I T I ON S , X < 0 > 



0 0 



INITIAL TIME = 0 



FINAL TIME = 4 



TIME STEP = 

v •-/ •./ •./ •./ ■./ •./ 

.■% . N . N A VS .-\ .*% /S 



0 1 



: ¥ ¥ : 



: ¥ ¥ ¥ 



FREQUENCY OF OUTPUT 
£ 4 4 4 » £ X X 4 



Figure 3 GTRESP Output 
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Figure 3 (Cont.) 



GRAPHICAL TIME RESPONSE 

PROBLEM IDENTIFICATION -> EXAMPLE 



T 


Y ( T) 


U<T) 


X 1 (T) 


X2CT) 


8 


9 


25 


3 


9 


.1 


2.19944778 


19.2237137 


.111344751 


2.33333393 


.2 


3.34344475 


14.3947742 


.39242274 


3.454344 


.3 


5.3245394 


13.4744595 


.782794787 


4.24374241 


.4 


5.81339897 


7.39775741 


1.22749837 


4.5359999 


.5 


4.28493932 


5.3453327 


1.4839279 


4.59401112 


.599999993 


4.51352439 


3.37443244 


2.13399243 


4.37142174 


.4 99999996 


4.55142535 


2.2243938? 


2.55813341 


3.99323444 


.799999994 


4.44312831 


1.51321138 


2.93459524 


3.52552595 


.399999991 


4.28325128 


1.14451514 


3.24194244 


3.31323342 


.999999939 


4.84733255 


1.33493374 


3.53319394 


2.53333841 


1.89999999 


3.42337344 


-13.1932884 


3.4549387? 


-.9345453122 


1.19999993 


1.49843535 


-13.1375299 


3.5544344 


-1.35738124 


1.29999998 


.258794742 


-8.91332349 


3.33441444 


-3.34742137 


1.39999998 


-.748358824 


-5.43739354 


2.94349324 


-3.73255737 


1.49999993 


-1.4543345 


-2.79448431 


2.57239891 


-4.92729341 


1.59999998 


-1.34437253 


-.759574434 


2.14799317 


-4.33197549 


1.49999998 


-2.35893243 


.714871183 


1.77333415 


-3.33223478 


1.7999999? 


-2.39253433 


1.71454348 


1.49593844 


-3.49844248 


1.39999997 


-2.31383574 


2.3-3453974 


1.37418935 


-3.38499541 


1.9999999? 


-1.35192291 


2.45573142 


.73943324 


-2.4415531? 


2.89999994 


-1.44447353 


2.7443357 


.547919244 


-2.19459277 


2.19999994 


-1.41983913 


2.47275259 


.24998815 


-1.74997733 


2.29999994 


-1.18722972 


2.43587785 


.192873534 


-1.33313331 


2.39999994 


-.944037342 


2.22394833 


.3724495393 


-1.33448455 


2.49999994 


-.753194579 


1.93544372 


-.3148449982 


-.74212947! 


2.59999995 


-.574947448 


1.4314455 


-.9774397299 


-.497314733 


2.49999995 


-.414853731 


1.33489444 


-.117134954 


-.299748744 


2.79999995 


-.284448821 


1.35347972 


-.139922944 


-.145425754 


2.39999995 


-.174314194 


.313734237 


-.147445513 


-.0293594731 


2.99999995 


-.392342342? 


.595121514 


-.145984544 


.353924521 


3.39999994 


-.3274993774 


.413224971 


-.137593411 


.199391233 


3.19999994 


.3198515238 


.24422934? 


-.124743259 


.14379478 


3.29999994 


.3513344439 


. 145923933 


-.19949421? 


.149442341 


3.39999994 


.3713353237 


.3552349217 


-.3938872394 


.14489241 


3.49999993 


.381533273 


-.3115345112 


-.3748934485 


.158398742 


3.59999993 


.3347935445 


-.35793434 


-.9414148385 


.144497353 


3.49999993 


.832838323 


-.3377345728 


-.8477441939 


.139582217 


3.79999993 


.3772893214 


-.134315384 


-.3355433192 


.112852441 


3.89999993 


.3494859387 


-.113744432 


-.9251881781 


.3944741143 


3.99999992 


.0434847343 


-.109494415 


-.3144881498 


.3773949274 



Figure 3 
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Figure 3 (Cont.) 



PROBLEM IDENTIFICATION - EXAMPLE 

XX GRAPHICAL TIME RESPONSE SS 




ABSCISSA -> TIME AXIS 
ORDINATE -> RESPONSE MAGNITUDE 
TIC MARKS SHOW INTERVALS OF UNITY 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, 0 TO 4 
ORDINATE, -5 TO 5 

SYSTEM RESPONSE 

VARIABLE SYMBOL 

X2 + 

ERRORX 

INPUTT 



Figure 3 
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IV. TRANSFER FUNCTION PROGRAMS 



The two transfer function programs discussed in this 
section are modified versions of programs of the same names 
found in [Ref. 1] . 

The Frequency Response program (FRESP), discussed in 
part A of this section, determines and plots the frequency 
response of a transfer function over a specified range of 
frequencies. The output may take the form of rectangular 
Bode plots or a polar Nyquist plot or both as desired. 

The Root Locus program (RTLOC), discussed in part B, 
calculates and plots the root locus of a transfer function 
for a specified range of gains. It is also possible to 
enlarge a small rectangular section of the root locus for 
more detail. 
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A. FRESP 



1 . Frequency Response Program (FRESP) 

a. Introduction 

The response of a system as a function of fre- 
quency is a very important characteristic of that system. 

Some common ways of graphically displaying the frequency 
response of a system include the use of amplitude and phase 
Bode plots and Nyquist diagrams. These are valuable tools 
for system analysis. 

b. Description of Program 

The FRESP program [Ref. 1: pp. 105-113] is 

used to determine the frequency response of a rational trans- 
fer function G(s) of the form 



G(s) 



A(s) 

B(s) 



where 

A(s) = a^ + agS + ... + a^s + s 

N-l N 

B(s) = b^ + b 2 s + ... + b^s + s 

The output may take the form of Bode plots or a Nyquist dia- 
gram or both in addition to tabular data. 

FRESP gives the user the option of supplying 
discrete frequency values or allowing the program to linearly 
or logarithmically interpolate frequency values between two 
limit values. The complex number G(jCJ) can be computed for 
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each frequency value. G(jGJ) can be written in rectangular 
form as 

G(jCJ) = R(CO) + jX(CJ) 

where R(CJ) and X(GJ) are real values. The magnitude and 
phase of G( j QJ) can then be written as 

G(jCJ) = [R 2 (6 l)> + X 2 (OJ)] 4 
arg G(jCJ) = arctan X(CJ)/R(CJ) 

The FRESP program uses the subroutines PROOT , 

PVAL, and SEMBL. The subroutines MAXI, PHNOM, GRAPH, and 
SPLIT that appear in the FORTRAN version do not appear in the 
BASIC version. They are incorporated into the main program 
and into various plotting subroutines written for the specific 
microcomputer system used. 

PROOT. This subroutine uses a modified Bairstow 
method for root extraction to determine the roots of a 
polynomial with real coefficients. It is used to determine 
the roots of the numerator and denominator of G(s) when 
they are entered in polynomial form. 

SEMBL. This subroutine determines the coeffi- 
cients of a polynomial from the roots of the polynomial. 

It is used to determine the polynomial coefficients of the 
numerator and denominator of G(s) when they are entered in 
the factored form. This subroutine together with PROOT 
provides the feature that G(s) may be entered in either 
factored or coefficient form and it will appear in the output 
in both forms. 
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PVAL. This subroutine is used to evaluate a 



polynomial A(s) with real coefficients for 
s = PR + jPI 

PVAL is executed twice for each frequency value used, once 
for the numerator of the transfer function and again for 
the denominator. 

c. Program Translation Problems 

In addition to the programing considerations dis- 
cussed in section II, an output anomaly was traced to an 
apparent oversight in the subroutine PVAL. It was discovered 
that a zero order numerator over a second order denomina- 
tor with a free s such as 
1/s ( s+10 ) 
was plotted as 

1/ (s+10) 

The reason for this was that PVAL treated a zero order 

polynomial and a first order polynomial both as a first order 

polynomial which in this case resulted in the free s in 

the denominator being cancelled. The reason for this can 

be seen by examining the portion of FORTRAN code from PVAL 

in [Ref. 1: p. 164] repeated below: 

P=CMPLX( A(NN+1 ) ,0. ) 

DO 100 J=1 ,NN 
100 P=P*S+A(NN+1-J ) 

VR=REAL(P) 

In this portion of code NN represents the order 
of the polynomial being evaluated and the intention is that 
the DO loop be executed a number of times equal to the order 
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of the polynomial. But because of the nature of the DO loop, 
if the polynomial is of zero order (i.e. NN equals zero) 
the loop will be executed once just as if the polynomial 
was of first order [Ref. 3]. One way to correct this 
problem is to modify the section of PVAL as shown here: 

P=CMPLX(A(NN+1) ,0. ) 

IF (NN) 110,110,90 
90 DO 100 J=1 ,NN 
100 P=P*S+A(NN+1-J ) 

110 VR=REAL(P) 

The essence of the modification is to skip the DO loop if 
the polynomial has order zero. Translating this modifica- 
tion into the BASIC version of FRESP solved the problem. 

2 . FRESP User 1 s Guide 

Push "return" after each input. 



STEP 1: 

SCREEN PROMPT: 

FREQUENCY RESPONSE PROGRAM (FRESP) 

THIS PROGRAM OBTAINS AND PLOTS THE FREQUENCY 
OF A RATIONAL TRANSFER FUNCTION OVER A SPECIFIED 
RANGE OF FREQUENCIES. BOTH RECTANGULAR BODE 
PLOTS AS WELL AS A POLAR NYQUIST PLOT CAN BE 
OBTAINED. 

PROBLEM IDENTIFICATION -> 



REMARKS : 

Label the problem. Type in any appropriate combina- 
tion of letters, numbers, and/or symbols, excluding 
commas and colons. This input is limited to 255 
characters . 
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STEP 2: 



SCREEN PROMPT: 

INPUT THE GAIN -> 

REMARKS : 

Enter the transfer function gain. 



STEP 3: 




SCREEN 


PROMPT : 




'KEY' = P FOR POLYNOMIAL COEFFICIENT FORM 
'KEY' = F FOR POLYNOMIAL FACTORED FORM 

INPUT 'KEY', ORDER FOR NUMERATOR ,A(S) -> 



REMARKS : 

Choose the preferred method of entering the numera 
tor polynomial of the transfer function, G(s). 

STEP 4A (assumes coefficient form chosen): 



SCREEN 


PROMPT : 

INPUT POLYNOMIAL COEFFICIENTS 




INPUT COEFFICIENTS IN ASCENDING PWRS OF S 
INPUT COEFFICIENT OF S~0 -> 



REMARKS : 

Enter the coefficient along with its algebraic sign 
as prompted. The program assumes that the coeffi- 
cient of the highest power of s is one. 

STEP 4B (assumes factored form chosen): 



SCREEN 


PROMPT : 

INPUT POLYNOMIAL FACTORS 

FOR FACTOR 1 
REAL PART -> 
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IMAGINARY PART -> 



REMARKS : 

Enter the factor with its algebraic sign. Do not 
enter the associated root. (e.g. For "(s-1)" 
enter "-1".) After the real part of the factor is 
entered, the program asks for the imaginary part 
allowing for input of quadradic factors. If a 
non-zero imaginary part is entered, the program 
automatically enters its complex conjugate. 

STEP 5: 

SCREEN PROMPT: 

•KEY' = P FOR POLYNOMIAL COEFFICIENT FORM 
'KEY' = F FOR POLYNOMIAL FACTORED FORM 

INPUT 'KEY', ORDER FOR DENOMINATOR, B(S) -> 

REMARKS : 

The denominator is entered just as described for 
the numerator in steps 3 and 4. 

STEP 6: 

SCREEN PROMPT: 

MINIMUM FREQUENCY VALUE -> 

REMARKS : 

Enter the minimum frequency value in radians per 
second of the frequency range of interest. 

STEP 7: 

SCREEN PROMPT: 

MAXIMUM FREQUENCY VALUE -> 
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REMARKS : 



Enter the maximum frequency value in radians per 
second of the frequency range of interest. 

STEP 8: 

SCREEN PROMPT: 

NUMBER OF FREQUENCY VALUES TO BE USED -> 

REMARKS : 

Enter an integer from 1 to 200. Since the horizon- 
tal resolution capability of the graphics on the 
microcomputer system used is 191 pixels, any number 
of values greater than 191 adds no more detail to the 
Bode plots. For the Nyquist diagram the greater the 
number of values the greater the detail due to the 
nature of the plot. If a list of discrete fre- 
quency values will be supplied (i.e. option 1 in 
step 9) enter the number of frequencies in that 
list . 

STEP 9: 

SCREEN PROMPT: 

0 = LOGARITHMIC INTERPOLATION 

1 = DISCRETE VALUES SUPPLIED 

2 = LINEAR INTERPOLATION 

CHOOSE ONE -> 

REMARKS : 

Make a choice by entering the associated number. This 
step allows the user to choose the method used by the 
program to obtain discrete frequency values used to 
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generate data points. LOGARITHMIC INTERPOLATION pro- 
duces data points that appear equally spaced on a 
logarithmic scale between the minimum and maximum 
frequency values entered in steps 6 and 7. With this 
choice a Bode and/or Nyquist plot may be obtained. 
LINEAR INTERPOLATION produces data points that appear 
equally spaces on a linear scale. Only a Nyquist 
plot may be obtained with this choice. DISCRETE VALUES 
SUPPLIED allows the user to supply a list of 
discrete frequency values of interest. No plots 
may be obtained, however, with this option. Below 
is a table summarizing the available graphical output 
for each choice. 



CHOICE 


BODE 


NYQUIST 


0 


YES 


YES 


1 


NO 


NO 


2 


NO 


'YES 



STEP 10 (assumes choice 1 in step 9): 

SCREEN PROMPT: 

TYPE IN THE DISCRETE FREQUENCY VALUES, 

(number entered in step 8) VALUES NEEDED 
FREQ (1)? 

REMARKS : 

Enter the list of frequency values as prompted. When 
the last value has been entered the program will run 
to completion outputting data similar to Figure 4 with 
the exception of the plots. 
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STEP 11 (assumes choice 0 in step 9): 

SCREEN PROMPT: 

BODE PLOT? (Y/N) -> 

REMARKS : 

The program outputs tabular data regardless of the 
choice made here. 

STEP 12 (assumes choice 0 or 2 in step 9): 

SCREEN PROMPT: 

NYQUIST PLOT? (Y/N) -> 

REMARKS : 

The program outputs tabular data regardless of the 
choice made here. 

After this step, the program runs to completion. 

See sample output in Figure 4. 

3 . FRESP Example 

It is desired to know the frequency response of the 
transfer function 

8(0.5 + s) 

G(s) = 

4 + 6s + 3s 2 + s 3 

The frequency range of interest is 6J = 0 . 1 rad/sec to 100 
rad/sec. A Bode plot of the magnitude and phase is desired 
but a Nyquist plot is not. The plots should be generated 
from 100 values of frequency logarithmically spaced from 
0.1 to 100. 
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Referring to the steps described in the previous 
section, FRESP User's Guide, this problem would be entered 
as follows: 



return . 



STEP 


1 . 


Enter "EXAMPLE CR " where <CR> 


STEP 


2. 


Enter "8 <CR>." 


STEP 


3. 


Enter "F,l <CR>." 


STEP 


4B . 


Enter ".5 <CR> 0 <CR>." 


STEP 


5. 


Enter "P , 3 <CR> 4 <CR> 6 <CR> 3 


STEP 


6. 


Enter ".1 <CR>." 


STEP 


7. 


Enter "100 <CR>." 


STEP 


8. 


Enter "100 <CR>." 


STEP 


9. 


Enter "0 <CR>." 


STEP 


11. 


Enter "Y <CR>." 


STEP 


12. 


Enter "N <CR>." 



The program will then continue to completion pro- 
ducing the output seen in Figure 4. 

The output begins with the program name and problem 
identification. The transfer function gain is given followed 
by the transfer function numerator and denominator each 
listed as coefficients in ascending powers of s and then as 
the real and imaginary parts of the roots. 

The second page of output is tabular data. It 
consists of the radian frequency and the transfer function's 
corresponding real and imaginary parts, magnitude, and phase 
in radians and in degrees. Although 100 frequency values 
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were generated and their associated data points plotted, 
tabular data appears only for every other data point. To 
reduce unnecessary output and for formatting purposes the 
following scheme is used for tabular data. 



NUMBER OF FREQUENCIES TABULAR DATA 

REQUIRED FROM STEP 8 PRINTED FOR 



0 to 50 
51 to 100 
101 to 150 
151 to 200 



every freq. 
every other freq. 
every 3rd freq. 
every 4th freq. 



The next two pages of output are the Bode plots for 
amplitude and phase. The plots are headed with the problem 
identification and type of plot and below each plot is 
sufficient information to allow proper interpretation of the 
plotted data. Note that the phase angles are normalized to 
always remain between -180 and +180 degrees. 

This is the end of the output generated from the 
example as it was input since a Nyquist plot was not desired. 
A Nyquist plot is included as the last page of Figure 4 for 
the sake of completeness. 
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FREQUENCY RESPONSE 
PROBLEM IDENTIFICATION - EXAMPLE 
****** * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * ****** * * * * * 



GAIN=3 



NUMERATOR COEFFICIENTS - IN ASCENDING POWERS OF S 

.5 

1 



NUMERATOR ROOTS ARE 
REAL PART 



IMAGINARY PART 



-.5 



0 



DENOMINATOR COEFFICIENTS - IN ASCENDING POWERS OF 



4 

6 

3 

1 



DENOMINATOR ROOTS ARE 
REAL PART 



IMAGINARY PART 



-1 

-1 

-1 



-1 .73285081 
1 .73205081 
0 



* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 



Figure 4 



FRESP Output 



Figure 4 (Cont.) PROBLEM IDENTIFICATION - EXAMPLE 
RADI/W FREQ. REAL PART IMAGINARY PART WGNITUDE PMSE(RAD) PHASE'! DEG) 



. 107224722 1.01795827 

. 123284474 1.92247107 

. 141747416 1.92954848 

.142975084 1.33385084 

.187381743 1.05094983 

.215443449 1.34445059 

.247707434 1.08493827 

.284893588 1.11285044 

.327454917 1.14575147 

.374493582 1.18495955 

.43287413 1.23774441 

.497702359 1.29903747 

.572234749 1.37100815 

.457933228 1.45241093 

.754443332 1.53955702 

.849749008 1.42457894 

1.80000001 1.49230749 

1.149757 1.71474725 

1.32194114 1.44344372 

1.51991109 1.4103491 

1.74752841 .945794053 

2.00923302 .330509559 

2.31012972 -.134534879 

2.45408781 -.423312809 

3.05335554 -.59737737 

3.51119177 -.473004414 

4.0370173 -.404413557 

4.44158888 -.32894124 

5.33449929 
4.13590734 
7.85488239 
3.1113084 
9.32403353 
10.7224724 
12.3284474 
14.1747413 
14.2975084 
18.7381745 
21.5443472 
24.7797439 
28.4883591 
32.7454921 
37.4493587 
43.2874135 
49.7702345 
57.2234774 
45.7933234 
75.4443342 
84.9749019 
100.000092 



.0514283748 1.01334782 

.0584342913 1.92415089 

.0443192482 1.93170223 

.074499234 1.94152594 

.9832923195 1.95424531 

.9920537744 1.97041542 

.199388178 1.99152705 

.197134157 1.11799545 

.111131334 1. 1511234 

.110214475 1.19284571 

.191419543 1.24189277 

.0394138318 1.38153489 

.8422474713 1.37145954 

-.0298852383 1.45254188 

-.118040545 1.54407711 

-.249401144 1.4453479 

-.441533444 1.75411404 

-.732787824 1.34478988 

-1.07434971 1.94392254 

-1.45914418 2.92239944 

-1.7451803 1.99459581 

-1.7944947 1.83454744 

-1.5559544 1.54174192 

-1.14877489 1.24324255 

-.393954307 .950479441 

-.52954775 .713394931 

-.342379243 .531334177 

.394282748 
.295723333 
.221178139 
.145791393 
.124515475 
.8934454855 
.8735491932 
.3531914145 
.049135442 
.3393024494 
.0228892471 
.8 172954174 
.0139724492 
9. 8324953 IE-93 
7.4721 1725E-83 
5.45830488E-03 
4.2738587E-83 
3.23173108E-03 
2.444294 15E-83 
1.84879722E-03 
1.3984 1954E-93 
1. 0577380 IE-03 
3.89129944E-94 
Figure 4 



.0597189222 


2.93593123 


.0572328479 


3.2320444 


.0443257541 


3.48559544 


.0714873348 


4.1373332 


.0790899987 


4.53144892 


.9840884104 


4.93250435 


.0929249474 


5.27275758 


.0959748558 


5.4990249 


.094491838 


5.54903421 


.9925904297 


5.3059542 


.3817543458 


4.48429433 


.3419889349 


3.55124725 


.0303193442 


1.74584839 


-.0143737124 


-.823339329 


-.3745349743 


-4.33513253 


-.159954477 


-9.11324809 


-.244252052 


-15.2551243 


-.403859573 


-23.1389419 


-.579143334 


-33. 1325894 


-.799399724 


-45.7979991 


-1.34534155 


-41.0395945 


-1.34289719 


-78.3424434 


-1.45704794 


-94.9413347 


-1.91344312 


-199.931333 


-2.13377179 


-122.254142 


-2.30504541 


-132.370573 


-2.44102701 


-139.348595 


-2.55902432 


-144.135483 


-2.43342774 


-151.132287 


-2.71143854 


-155.3455 


-2.77251189 


-153.353237 


-2.32375481 


-141.739291 


-2.34722805 


-144.279447 


-2.98430789 


-144.404593 


-2.9348947 


-148.224909 


-2.94344104 


-149.792724 


-2.98792444 


-171.143947 


-3.3974041 


-172.311739 


-3.3259474 


-173.322511 


-3.94033484 


-174.193418 


-3.95359432 


-174.953129 


-3.84519255 


-175.417503 


-3.97599595 


-174.190883 


-3.93377724 


-174.487435 


-3.99132094 


-177.119738 


-3.39787757 


-177.495374 


-3. 19357722 


-177.321939 


-3.10853252 


-178.105858 


-3.11284111 


-173.352722 


-3.11458744 


-173.547383 



-.259108135 

-.201047424 

-. 154424892 

-.118278908 

-.099141942 

-.0485723775 

-.3528722415 

-.9395004317 

-.029941495 

-.8224335034 

-.0171782883 

-.9130954894 

-9.84425437E-83 

-7.45924918E-03 

-5.437817 14E-93 

-4.24591907E-03 

-3.22744824E-93 

-2.44194098E-03 

-1. 34744 144E-93 

-1. 3974554 IE-03 

-1.95734283E-03 

-7.99879835E-04 



-.229992452 

-.142544144 

-.8921935998 

-.05981042 

-.9339127442 

-.0253788175 

-.0145334577 

-.3198533543 

-7.1 1247 154E-03 

-4.44519419E-83 

-3.0422238E-83 

-2.0111427SE-33 

- 1.32 1427 14E-03 

-8.48521474E-34 

-5.70934204E-04 

-3.75450085E-34 

-2.44911144E-34 

-1.42394323E-84 

-1.34818434E-04 

-7.02459939E-35 

-4.42235779E-35 

-3.84084457E-35 

-2.9905 1985E-05 
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Figure 4 (Cont.) 



PROBLEM IDENTIFICATION - EXAMPLE 

XX BODE PLOT (AMPLITUDE) XX 




ABSCISSA -> COMMON LOG OF FREQUENCY 
ORDINATE -> COMMON LOG OF AMPLITUDE 

MINIMUM FREQUENCY SHOWN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOWN ON ABSCISSA = 100 RADIANS/SEC 
AMPLITUDE LIMITS OF BODE PLOT ARE +-88 DECIBELS 



Figure 4 
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Figure 4 (Cont.) 



PROBLEM IDENTIFICATION - EXAMPLE 

XX BODE PLOT < PHASE) XX 




ABSCISSA -> COMMON LOG OF FREQUENCY 
ORDINATE -> PHASE < DEGREES) 

TIC MARKS SHOW MULTIPLES OF 90 DEGREES 
MINIMUM FREQUENCY SHOWN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOWN ON ABSCISSA = 100 RADIANS/SEC 
MAXIMUM PHASE ON ORDINATE SCALE = 90 DEGREES 
MINIMUM PHASE ON ORDINATE SCALE = -180 DEGREES 



Figure 4 



Figure 4 (Cont.) 



PROBLEM IDENTIFICATION - EXAMPLE 

XX NYQUIST PLOT XX 




ABSCISSA -> REAL PART OF G<JW> 
ORDINATE -> IMAGINARY PART OF G<JW> 
TIC MARKS SHOW INTERVALS OF UNITY 
AXES CROSS AT ORIGIN 



Figure 4 
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B. RTLOC 



1 . Root Locus Program (RTLOC) 

a. Introduction 

When analyzing a system it is helpful to know the 
location of the closed loop poles in the s-plane. These 
closed loop poles are the roots of the characteristic 
equation and determine the basic characteristics of the 
transient response of a closed loop system. A root locus 
plot is a plot of the roots of the characteristic equation, 
usually as a function of the gain of the transfer function, 
and therefore it is a valuable tool for system analysis. 

b. Description of Program 

The RTLOC program [Ref. 1: pp. 114-121] calcu- 

lates and plots the roots of the equation 
1 + KG( s ) = 0 

as a function of K. G(s) is assumed to be a rational function 
of the form 



G(s) 



N(s) 

D(s) 



and the root locus becomes the locus of roots of D^.(s) as K 
varies where 



D k (s) = D(s) + KN(s) 



RTLOC uses an algebraic plus linear progression 
scheme to vary K to give reasonable spacing of the roots. 

D^.(s) is obtained for each value of K and the subroutine PROOT 
is used to calculate its roots. 
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The scheme used to calculate values of K assumes 



that K takes on only positive values. If K is to range 
through negative value, the value of smaller magnitude (less 
negative) must be used as the minimum value. The routine 
starts at the maximum (less negative) value and becomes 
increasingly negative until the lower limit is reached. 

If both positive and negative values are desired then two 
separate runs must be made with only positive and only 
negative values. 

The RTLOC program has a feature that allows the 
user to specify a range of s and Cl) values around an area 
of the root locus plot that is of interest . That area is 
then enlarged and more closely spaced values of K are generated 
giving more detail to the plot. 

The major subroutines used are PROOT and SEMBL. 

The subroutine SPLIT used in the FORTRAN version is replaced 
by additional coding in the main program and by various 
plotting subroutines. 

PROOT. This subroutine uses a modified Bairstow 
method for root extraction to determine the roots of a 
polynomial with real coefficients. It is used to determine 
the roots of the numerator and denominator of G(s) when it 
is entered in polynomial form. It is also used to determine 

the roots of D-^s) for each value of K generated. 

K. 

SEMBL. This subroutine determines the coeffi- 
cients of a polynomial from the roots of the polynomial. It 
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is used to determine the polynomial coefficients of the 
numerator and denominator of G(s) when they are entered in 
the factored form. This subroutine together with PROOT 
provides the feature that G(s) may be entered in either 
factored or coefficient form. It appears in the output in 
both forms. 

c. Program Translation Problems 

In addition to the programing considerations 
discussed in section II, two problems were encountered. 

When a range of negative gains is entered the program cal- 
culates only one value of K due to the logic of line 229 in 
the FORTRAN code. This works correctly only for positive 
ranges of gain. (Note that K is represented by the variable 
G in the computer program. ) In the BASIC version coding 
was added to test the sign of the range of gains and modify 
the logic to correctly handle the negative case. 

The second problem was encountered when using 
gains of very small magnitude necessary in using this 
program for w' -plane analysis. The schemes used to generate 
values of K for a reasonable spacing of the roots when plotted 
worked fine for the magnitudes usually encountered in s-plane 
analysis but was not flexible enough to handle gains of much 
smaller magnitudes (e.g. 0 to IE-3). To add this necessary 
flexibility the following FORTRAN lines used to generate 
values of gain G, 
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227 G=1 . 15* (G+SIGNG*0 . 05 ) 

228 G=1.04*(G+SIGNG*0.02) 

were modified by replacing the constant values 0.05 and 0.02 

with the variables D1 and D2 , respectively, where 

D1=ABS ( GMIN-GMAX ) / 700 
D2=ABS(GMIN-GMAX) /1500 

so that the gain increment was a function of the range of 
gains of interest . 

2 . RTLQC User's Guide 

Push "return" after each input. 

STEP 1: 

SCREEN PROMPT: 

ROOT LOCUS PROGRAM ( RTLOC ) 

THIS PROGRAM PLOTS THE ROOT LOCUS OF A 
DESCRIBED SYSTEM. 

PROBLEM IDENTIFICATION -> 

REMARKS : 

Label the problem. Type in any appropriate combina- 
tion of letters, numbers, and/or symbols, excluding 
commas and colons. This input is limited to 255 
characters . 

STEP 2: 

SCREEN PROMPT: 

INPUT THE RANGE OF GAINS (MIN, MAX) -> 

REMARKS : 

Enter the range of transfer function gains of in- 
terest. Note that the minimum gain is the gain of 
lowest absolute value. 
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STEP 3: 



SCREEN 


PROMPT : 

'KEY' = P FOR POLYNOMIAL COEFFICIENT FORM 
'KEY' = F FOR POLYNOMIAL FACTORED FORM 

INPUT 'KEY', ORDER FOR NUMERATOR, N(S) “> 



REMARKS : 

Choose the preferred method of entering the 
numerator polynomial of the transfer function 
G(s) . 

STEP 4A (assumes coefficient form chosen): 



SCREEN 


PROMPT : 

INPUT POLYNOMIAL COEFFICIENTS 

INPUT COEFFICIENTS IN ASCENDING PWRS OF s OR 




w * 

INPUT COEFFICIENT OF S~0 -> 



REMARKS : 

Enter the coefficient along with its algebraic sign 
as prompted. The program assumes that the 
coefficient of the highest power of s is one. 

STEP 4B (assumes factored form chosen): 



SCREEN 


PROMPT : 

INPUT POLYNOMIAL FACTORS 

FOR FACTOR 1 
REAL PART -> 




IMAGINARY PART -> 
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REMARKS : 



Enter the factor with its algebraic sign. Do not 
enter the associated root. (e.g. For "(s-1)" enter 
After the real part of the factor is 
entered the program asks for the imaginary part 
allowing for input of quadradic factors. If a 
non-zero imaginary part is entered the program 
automatically enters its complex conjugate. 

STEP 5: 

SCREEN PROMPT: 

’KEY' = P FOR POLYNOMIAL COEFFICIENT FORM 
'KEY' = F FOR POLYNOMIAL FACTORED FORM 

INPUT 'KEY', ORDER FOR DENOMINATOR, D(S) “> 

REMARKS : 

The denominator is entered just as described for the 
numerator in steps 3 and 4. 

STEP 6: 

SCREEN PROMPT: 

WOULD YOU LIKE TO LOOK AT ONLY A PART OF THE 
ROOT LOCUS? (Y/N) -> 

REMARKS : 

Enter an "N" for the option of viewing the root locus 
for the entire range of gains entered in STEP 2. 

Enter a "Y" for the option of viewing only a portion 
of the root locus defined by minimum and maximum 
values of sigma and CJ of interest to the user. 
Usually, it is best to view the entire root locus 
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first then decide if any section needs to be 
enlarged to reveal more detail and rerun RTLOC 
to view only that portion. 

STEP 7 (assumes "Y" entered in step 6): 

SCREEN PROMPT: 

ENTER SIGMA MIN, SIGMA MAX -> 

REMARKS : 

Enter the minimum and maximum values of sigma (i.e. 
the real axis) of the portion of the root locus 
plot of interest. 

STEP 8 (assumes "Y" entered in step 6): 

SCREEN PROMPT: 

ENTER OMEGA MIN, OMEGA MAX -> 

REMARKS : 

Enter the minimum and maximum values of omega (i.e. 
the imaginary axis) of the portion of the root 
locus plot of interest. 

STEP 9: 

SCREEN PROMPT: 

PRINT OUT OF GAIN DATA? (Y/N) -> 

REMARKS : 

This gives the user the option to suppress the print 
out of the gain data if it is not needed. This saves 
time and paper if only the root locus plot is desired. 
After this step the program runs to completion. 
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3 . RTLOC Example 



It is desired to know the root locus for the open 
loop transfer function G(s) given by 



G(s ) 



1 . 2 + s 

s(8 + 9s + s^) 



for gain K of from 0 to 30. Referring to the steps des- 
cribed in section 2, RTLOC User's Guide, this problem would 
be entered as follows: 

STEP 1. Enter "EXAMPLE <CR>" where <CR> = carriage 

return . 



STEP 2. Enter "0,30 <CR>." 

STEP 3. Enter "F,l <CR> . " 

STEP 4B . Enter "1.2 <CR> 0 <CR>." 

STEP 5. Enter "F,3 <CR> 0 <CR> 0 <CR> 1 <CR> 0 <CR> 
8 <CR> 0 <CR> . " 



STEP 6. Enter "N <CR> . " 

STEP 9. Enter "Y <CR>." 

The program will then run to completion producing 
the output seen in Figure 5. 

The output begins with the program name and problem 
identification. Next the transfer function numerator and 
denominator are listed as coefficients in ascending powers 
of s and then as the real and imaginary parts of the roots 
(i.e. open loop zeros for the numerator and open loop poles 
for the denominator). The minimum and maximum values of gain 
are listed next. 
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If the option to view only a part of the root locus 
is selected, this fact is stated next followed by the ranges 
of sigma and omega that are chosen. (See Figure 6) 

The next several pages of output contain the gain 
data which may be suppressed with the appropriate response in 
step 9. The data point number is listed with the value of 
gain and the real and imaginary parts of the corresponding 
roots of the open loop system. 

The last page of output contains the root locus plot 
itself. It begins with the problem identification and 
heading and ends with sufficient data listed to allow proper 
interpretation of the plot. 

Figure 6 shows the output from the same example 
problem with the gain data suppressed and the option to 
view only a part of the root locus chosen. 
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ROOT LOCUS 

PROBLEM IDENTIFICATION - 

X X X X X X X X X X X X X X X X X X X X X X X X 



EXAMPLE 

X X X X X X X X X 



X X X X X X X X X X X X X X X X X 



NUMERATOR COEFFICIENTS - IN ASCENDING POWERS OF S OR W' 



1.2 

1 



OPEN-LOOP ZEROS 

REAL PART IMAGINARY PART 

- 1.2 0 



DENOMINATOR COEFFICIENTS - IN ASCENDING POWERS OF S OR W' 

0 

8 

9 

1 



OPEN-LOOP POLES 
REAL PART 

-1 

-8 

0 



IMAGINARY PART 

0 

0 

0 



MIN GAIN MAX GAIN 

0 30 



xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 



Figure 5 RTLOC Output 



76 



Figure 
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PART 








-1 


0 










-3 


0 










0 


0 




GAIN 


= 


.9492357143 


ROOTS 


ARE 










REAL PART 


I MAG. 


PART 








-7.4039 1563E-03 9 










-.993530 1 1 1 


0 










-7.9940 1093 


0 




GAIN 


= 


. 105964236 


ROOTS 


ARE 










REAL PART 


I MAG. 


PART 








-.9 159695125 


0 










-.996917641 


0 










-7.93711234 


0 




GAIN 


= 


. 17 1 144643 


ROOTS 


ARE 










REAL PART 


I MAG. 


PART 








-.0258639696 


0 










-.994965313 


0 










-7.97916572 


0 




GAIN 


= 


.246102054 


ROOTS 


ARE 










REAL PART 


I MAG. 


PART 








-.037328049 


0 










-.99 2664433 


0 










-7.97000752 


9 




GAIN 


= 
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Figure 5 
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Figure 5 (Cont.) 
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Figure 5 
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Figure 
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Figure 5 
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Figure 5 (Cont.) 



PROBLEM IDENTIFICATION - EXAMPLE 

XX ROOT LOCUS PLOT XX 




ABSCISSA -> REAL (SIGMA) AXIS 
ORDINATE -> IMAGINARY (OMEGA) AXIS 
TIC MARKS SHOW INTERVALS OF 1 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -8 TO 1 
ORDINATE, -4 TO 5 



Figure 5 
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ROOT LOCUS 

PROBLEM IDENTIFICATION - EXAMPLE WITH OPTION 

X * X * X X * * * * * ****** * X * ***** X * * * * * * * X X £ X X X 



NUMERATOR COEFFICIENTS - IN ASCENDING POWERS OF S 



1.2 

1 



OPEN-LOOP ZEROS 
REAL PART 



IMAGINARY PART 



- 1.2 



0 



DENOMINATOR COEFFICIENTS - IN ASCENDING POWERS OF 

0 

S 

? 

1 



OPEN-LOOP POLES 
REAL PART 

-1 

-3 

0 



IMAGINARY PART 

0 

0 

0 



MIN GAIN MAX GAIN 

0 30 



* X * X x X X XXXXXXXXXXXX x X X X X X X X X X X X X X X X X X X X X X X X X X X X X X 



OPTION TAKEN 

SIGMA MIN = -2 SIGMA MAX = 0 

OMEGA MIN = -1 OMEGA MAX = 1 

Figure 6 RTLOC Output with Option 
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OR W' 



S OR W' 



Figure 6 (Cont.) 



PROBLEM IDENTIFICATION - EXAMPLE WITH OPTION 
XX ROOT LOCUS PLOT XX 




ABSCISSA -> REAL < SIGMA) AXIS 
ORDINATE -> IMAGINARY < OMEGA) AXIS 
TIC MARKS SHOW INTERVALS OF 1 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -2 TO 1 
ORDINATE, -1 TO 2 



Figure 6 
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V. w' -PLANE ANALYSIS 



Digital control systems are becoming more and more common 
Digital control laws have unique characteristics that can 
only be approximated by using classical techniques in the 
continuous s domain. 

In the w' domain all analog control system design tech- 
nology transfers completely for digital control system design 
An important advantage of the w' domain is that non-minimum 
phase effects of the sampling and data-hold operations and 
of sampling rate can be directly accounted for without 
approximation while using conventional frequency domain 
design and analysis tools such as root locus and Bode 
plots. [Ref. 4] 
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A . BACKGROUND 



This s domain is used for continuous system analysis. 
When a digital system is considered the z, w, or w' domain 
must be used. These domains are related as shown below. 



z-1 sT 

w = = tanh — 

z+1 2 



w' = (2/T) w 



where T is the sampling period and 



tanh x = 



x x 
e - e 



x , x 
e + e 



In an s-plane root locus plot the region of stability is 
the left half plane, that is all roots with negative real 
parts. This region of stability is mapped into a unit circle 
with its center at the origin in the z-plane. By use of 
the bilateral transformation shown above, the stability region 
of the z-plane is mapped back into the left half plane to 
form the w-plane. The w' -plane takes it a step farther by 
multiplying the w-plane by the factor 2/T where T is the 
sampling period. This gives the w'-plane the property 
that not only is the region of stability the left half plane 
as in the continuous s-plane but w' approaches s as the sam- 
pling period T approaches zero [Ref. 5]. 
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B. TRANSFER FUNCTION APPROACH USING THE w' -PLANE 



In this section the second order system 



4* ~ ... 

^ V V 


10 






v ^ 


s ( s+10 ) 




r 



is analyzed using the two transfer function programs dis- 
cussed in part IV of this thesis. It is then converted 
to the w'-plane and analyzed for periods of .001 seconds, 

.01 seconds, and .1 seconds using the same two transfer func- 
tion programs. These results are used to gain insight into 
w'-plane analysis. 

The system is converted to the w'-plane by adding a 
sampler and a digital to analog converter in the form of a 
zero order hold. This modified system is shown below. 




The new open loop transfer function is transformed to 
the z-plane and then to the w'-plane as follows: 



G(s) = 10K( l-e** sT ) 

s^ (s+10 ) 
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G(z) = 10K( 1-z 1 ) — 

10 



Tz 



-10T . 

( 1-e )z 



(z-l)‘ 



-10T 

10(z-l) (z-e XUi ) 



= K 



1-e 



-10T 



1n/ -10T , 

z-1 10(z-e ) 



For the period 



T= . 001 



G(z) reduces to 



G(z) = 



4.98337 x 10 6 K(z+. 996681) 
(z-l)(z-. 990049834) 



( 1 ) 



Similarly, for 
T= . 01 

4.837418 x 10" 4 K( z+ . 9672185 ) 

G(z) = 

(z-l)(z-. 904837418) 



And for 



T= . 1 

.03678794412 K( z+ . 718281827 ) 

G( Z ) = 

(z-l)(z-. 3678794412) 

Now G(z) is converted to the w' domain using the rela- 
tionship 
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1 + T/2 w' 

z = 

1 - T/2 w* 

For period 

T=. 001 

1 + .0005 w' 

z = 

1 - .0005 w' 



Substituting into equation (1) gives 



4.98337 x 10 6 K 



1+ . 0005w 1 



1- . 0005w 1 



G(w' ) = 



+ . 996681y 



1+.0005W' 



1-.0005W' 



- 1 



l+.0005w' 



1- . 0005w 1 



- .990049834 



which reduces to 



-4.156878 x 10 9 K (w' + 1202819) (w 1 - 2000) 

GCw') = (2) 

w' (w 1 + 9.9999) 

Note that the gain is 1.0K which is the same as the s domain 
gain. Similarly, for 

T= . 01 



l+.005w' 

z= 

1-.005W' 
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G(w ' ) 



(3) 



-4.1625027 x 10 6 K(w ' +12002 . 00418 ) (w ' -200 ) 
w' (w' + 9.991674985) 



And for 



T= . 1 



l+.05w' 

z = 

1- . 05w ' 

-.0037882843 K(w' + 121 . 9858708 ) (w ' -20 ) 

G(w') = (4) 

w' (w' + 9.242343139) 

Note again that the gain is 1.0K also in equations (3) and 
(4). 

1 . Frequency Response 
a. s-Plane 

The open loop transfer function of the system 
in the s domain is 

10K 

G(s) = 

s(s+10 ) 

The frequency response of G(s) with a gain K of 5 is described 
by the output of the FRESP program shown in Figure 7. From 
the Bode plots it is seen that the gain margin is infinite 
and the phase margin is 65.1 degrees. 
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b. w' -Plane with a Period of .001 Seconds 



The open loop transfer function for the equiva- 
lent sampled system is represented by equation (2) for a 
period T of .001 seconds. This transfer function is 
entered into the FRESP program as if w' were an s. The 
resulting output in Figure 8 is interpreted as described 
earlier in the thesis keeping in mind that all references 
to s are actually referring to w'. From this output it 
is seen that the gain margin is no longer infinite as it 
was in the continuous case but it is still quite high at 50 
dB. The phase margin has dropped slightly from 65.1 degrees 
to 64.6 degrees. 

c. w' -Plane with a Period of .01 Seconds 

If the sampling rate of the system is decreased 
so that the sampling period is .01 seconds, then the open 
loop transfer function is represented by equation (3). 
Inputting the transfer function into FRESP results in the 
output shown in Figure 9. From the Bode plots the gain 
margin is found to be 32 dB. This is down from 50 dB for 
the case of a period of .001 seconds indicating a decrease 
in stability. The phase margin is 64.5 degrees which is also 
lower than the previous case but only by one tenth of a 
degree . 

d. w'-Plane with a Period of .1 Seconds 

If the period is again increased to .1 seconds, 
the open loop transfer function is represented by equation 
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(4). Entering this transfer function into FRESP produces 
the output seen in Figure 10. From this output the gain 
margin is found to have decreased to 13.2 dB and the phase 
margin to 53.7 degrees. 

e. Summary of Frequency Response Results 

The table below is a brief summary of the gain 
and phase margins found in each case. 

CASE GAIN MARGIN PHASE MARGIN 



s-plane 
w' , T= . 001 
w' , T= . 01 
w' , T= . 1 



infinite 
50 dB 
32 dB 
13.2 dB 



65.1 deg. 

64.6 deg. 
64.5 deg. 

53.7 deg. 



It can be seen from this table that the con- 
tinuous case is the most stable. The sampled cases become 
less stable as the sampling period increases. It is also 
noticed that the gain margin is more sensitive to changes 
in the sampling period than the phase margin. 

2 . Root Locus 

Now a comparison is made using the root locus program 
using the same cases used above, 
a. s-Plane 

If the open loop transfer function 



10K 

G(s) = 

s(s+10 ) 



is entered into the RTLOC program, the result is the output 
in Figure 11. The root locus plot shows that the system 
never becomes unstable at any gain. 
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b. w' -Plane with a Period of .001 Seconds 



If the equivalent sampled transfer function with a 
sampling period of .001 seconds represented by equation (2) 
is entered into RTLOC , the output will be that found in 
Figure 12. From the root locus plot it can be seen that 
there is a slight tendency for the plot to curve toward the 
right half plane as it moves further from the real axis. The 
tendency is so slight that it will take a very large gain 
to drive the system unstable. 

c. w' -Plane with a Period of .01 Seconds 

If the sampling period is increased to .01 seconds 
the transfer function is represented by equation (3). 

Entering this transfer function into RTLOC results in the 
output seen in Figure 13. It can be seen that the tendency 
for the curve to bend toward the unstable right half plane 
is increased indicating that a lowe'r value of gain than in 
the previous case will drive the system to instability. 

d. w' -Plane with a Period of .1 Seconds 

When the sampling period is further increased to .1 
seconds represented by equation (4), the resulting RTLOC 
output is that seen in Figure 14. In the root locus plot 
for this case the tendency to become unstable is much more 
pronounced. Here the plot bends into the right half plane 
within the limited boundaries of the portion plotted. It 
becomes unstable for values of gain K greater then 24.6. 
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e. Summary of Root Locus Comparison 

It can be seen that the continuous system 

10K 

G(s) = 

s(s+10) 

is always stable. When this system is sampled it can be 
seen by equations (2), (3), and (4) that a zero is added 
in the right half plane making it a non-minimum phase 
system. It can also be seen that by decreasing or increasing 
the sampling period the distance of this zero from the 
origin increases or decreases respectively. The closer 
this zero is to the imaginary axis, the greater effect it 
has on bending the root locus into the right half plane. 

This effect can be seen by examining the root locus plots 
in Figures 11 through 14. 
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FREQUENCY RESPONSE 

PROBLEM IDENTIFICATION - EX 1 S-PLANE 
************************************************** 

GAIN=50 



NUMERATOR COEFFICIENTS - 


IN ASCENDING POWERS OF S 


1 





DENOMINATOR COEFFICIENTS 
0 

10 

1 


- IN ASCENDING POWERS OF S 


DENOMINATOR ROOTS ARE 
REAL PART 


IMAGINARY PART 


0 

-10 


0 

0 



************************************************** 
Figure 7 s-Plane Frequency Response Example 
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Figure 7 (Cont.) 

6 'PROBLEM IDENTIFICATION - EX 1 S-PLANE 



RADIAN FREQ. 


REAL PART 


IMAGINARY PART 


MAGNITUDE 


PHASE (RAD) 


PHASE(DEG) 


. 187224722 


-.499942519 


-44.4248044 


44.4274849 


-1.58151844 


-90.4143755 


. 123284674 


-.499924014 


-40.5503782 


40.5534598 


-1.58312422 


-90.7043485 


.141747414 


-.499899558 


-35.2449254 


35.2704484 


-1.58497017 


-93.3121334 


.142975884 


-.49984723 


-30.4713897 


30.4754427 


-1.58709244 


-90.9337311 


.187381743 


-.499824582 


-24.4741383 


24.4788128 


-1.58953234 


-91.973528 


.215443449 


-.499748927 


-23.1971749 


23.2025599 


-1.59233739 


-91.2342445 


.247707434 


-.499493392 


-20.1727384 


20.1788944 


-1.59554207 


-91.4190054 


.284803588 


-.499594743 


-17.54173 


17.5488429 


-1.59924984 


-91.4313988 


.327454917 


-.49944444 


-15.2529223 


15.2410977 


-1.48353017 


-91.8755439 


.374493582 


-.499292244 


-13.2414408 


13.2710345 


-1.48842794 


-92.1541447 


.43287413 


-.499844843 


-11.5290451 


11.5398417 


-1.41405498 


-92.4734859 


.497782359 


-.498744522 


-10.8213413 


10.3337455 


-1.42052558 


-92.3493393 


.572234749 


-.498348849 


-8.70912349 


8.72337108 


-1.42795771 


-93.2751394 


.457933228 


-.497844948 


-7.5448005 


7.58314024 


-1.43449501 


-93.7442907 


.754443332 


-.497155895 


-4.57289774 


4.5908749 


-1.44429891 


-94.324013 


.849749008 


-.49424488 


-5.78542399 


5.72714331 


-1.45755295 


-94.9798225 


1.00000001 


-.495049583 


-4.95049503 


4.97518594 


-1.47844503 


-95.7104301 


1.149757 


-.49347453 


-4.29200717 


4.32828294 


-1.43524942 


-94.5588594 


1.32194114 


-.491412428 


-3.71735478 


3.74949502 


-1.70222842 


-97.5305392 


1.51991109 


-.488710141 


-3.2153845 


3.25231425 


-1.72143305 


-98.4423427 


1.74752841 


-.485183205 


-2.77439444 


2.81847143 


-1.74380222 


-99.9125431 


2.00923382 


-.480598148 


-2.39194839 


2.43975239 


-1.74907957 


-101.340829 


2.31812972 


-.474448375 


-2.05472407 


2.10884075 


-1.79782443 


-103.897915 


2.45408781 


-.447953589 


-1.75841517 


1.81938454 


-1.83041085 


-104.374854 


3.85385554 


-.457347488 


-1.49740721 


1.54588448 


-1.84718449 


-104.981955 


3.51119177 


-.445123138 


-1.24772445 


1.34340148 


-1.90844789 


-109.347194 


4.0379173 


-.429931824 


-1.04497394 


1.14848199 


-1.95448982 


-111.984058 


4.44158888 


-.411372483 


-.884275139 


.977093108 


-2.03534225 


-114.398834 


5.33449929 


-.33914454 


-.729223288 


.324548597 


-2.34101572 


-118.087544 


4.13590734 


-.343241793 


-.591993409 


.494550959 


-2.12114917 


-121.532939 


7.05488239 


-.333844701 


-.473214233 


.579125104 


-2.18519089 


-125.29224 


8.1113984 


-.381580297 


-.371302282 


.479735431 


-2.25228744 


-129.344423 


9.32483358 


-.247415549 


-.284740944 


.392084033 


-2.32133517 


-133.902755 


10.7224724 


-.232584424 


-.21498901 


.318033049 


-2.39105392 


-134.997347 


12.3284474 


-.198419494 


-.140944331 


.2554847 


-2.44010143 


-140.953479 


14.1747418 


-.144155293 


-.117219273 


.203341927 


-2.5271982 


-144.797843 


14.2975984 


-.134758201 


-.0839135431 


.140450278 


-2.59123991 


-148.447144 


18.7381745 


-.110835455 


-.0591495479 


.125431075 


-2.45137337 


-151.912553 


21.5443472 


-.8884275129 


-.0411372448 


.8977093894 


-2.70702484 


-155.101248 


24.7787439 


-.0700481732 


-.028284442 


.0755424445 


-2.75739927 


-158.014845 


28.4893591 


-.0548748595 


-.0192483173 


.0581413881 


-2.8039212 


-140.452908 


32.7454921 


-.0424523898 


-.0138254234 


.0445949597 


-2.84520239 


-143.918147 


37.4493587 


-.0329494891 


-8.75147342E-03 


.0340919437 


-2.88197823 


-145.125248 


43.2874135 


-.0253314228 


-5.35193332E-03 


.0259987738 


-2.91454245 


-144.992187 


49.7792345 


-.0194018304 


-3.89827973E-03 


.0197895833 


-2.94338952 


-148.439273 


57.2234774 


-.0148147935 


-2. 5892747 IE-03 


.0150413339 


-2.94858484 


-170.887559 


45.7933234 


-.0112898382 


-1.71595499E-03 


.0114194984 


-2.99075403 


-171.357759 


75.4443342 


-8.58757143E-03 


-1.13522438E-03 


8.44228157E-03 


-3.01014044 


-172.449543 


84.9749019 


-4. 5234490 7E-03 


-7.58040404E-04 


4.54444572E-03 


-3.02711944 


-173.441243 


100.000002 


-4.95049495E-03 


-4.95049488E-04 


4.97518584E-03 


-3.04192405 


-174.239472 



Figure 7 
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Figure 7 (Cont.) 



PROBLEM IDENTIFICATION - EX 1 S-PLANE 

XX BODE PLOT (AMPLITUDE) XX 




ABSCISSA -> COMMON LOG OF FREQUENCY 
ORDINATE -> COMMON LOG OF AMPLITUDE 

MINIMUM FREQUENCY SHOWN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOWN ON ABSCISSA = 100 RADIANS/SEC 
AMPLITUDE LIMITS OF BODE PLOT ARE +-68 DECIBELS 



Figure 7 
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Figure 7 (Cont,) 



PROBLEM IDENTIFICATION - EX 1 S -PLANE 

** BODE PLOT (PHASE) XX 




ABSCISSA -> COMMON LOG OF FREQUENCY 
ORDINATE -> PHASE (DEGREES) 

TIC MARKS SHOW MULTIPLES OF 90 DEGREES 
MINIMUM FREQUENCY SHOWN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOWN ON ABSCISSA = 100 RADIANS/SEC 
MAXIMUM PHASE ON ORDINATE SCALE = -90 DEGREES 
MINIMUM PHASE ON ORDINATE SCALE = -180 DEGREES 



Figure 7 
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FREQ U EN CY R E S P ON 3 E 

PROBLEM IDENTIFICATION - EX 1 W' -PLANE T= .00 1 EXTENDED 

X X X X X X x x x x X x X X X x X X X XXXX X X X X X X X ¥ X X x X x x x x X X X X X X x X X X X 

GAIN=-2 . 973439E-08 

NUMERATOR COEFFICIENTS - IN ASCENDING POWERS OF S 

-2 . 495638E + 0’? 

120031? 

1 



NUMERATOR ROOTS ARE 
REAL PART 



IMAGINARY PART 



-120281? 8 

29 90 9 



DENOMINATOR COEFFICIENTS - IN ASCENDING POWERS OF S 



9 

?.???? 

1 



DENOMINATOR ROOTS ARE 
REAL PART 



IMAGINARY PART 



c o o 9 



a 

3 



X X X x X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X X XXXX X X 



Figure 8 w' -Plane Frequency Response Example, T= 



.001 
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Figure 



RADIAN FREQ. 

. 139749877 

.132194115 

.159228279 

.191791024 

.23181297 

.273255941 

.335140244 

.403701727 

.48424314 

.535782084 

.705439234 

.349753439 

1.02353103 

1.23234475 

1.43494327 

I. 73344954 
2.1544347 
2.59592423 
3.12571537 
3.74493533 
4.53437354 
5.44227724 
4.5793323 
7.72432905 
9.54543445 

II. 4975701 
13.3433433 
14.4810055 
20.9923392 
24.2012329 
29.1595339 
35.1119177 
42.2924292 
50.94 13S08 
41.3599735 
73.9072212 
39.9215997 
107.224724 
129.154943 
155.547414 
187.281745 
225.781975 
271.358323 
327.454921 
394.423412 
475.331923 
572.234775 
439.241221 
330.217532 
1000.30982 



< Cont ^PROBLEM 



I DENT I F I CAT I ON 



EX 1 M' -PLANE T=.0ul 



REAL PART 


IMAGINARY PART 


MAGNITUDE 


PHASE (RAD) 


PHASE (DEG) 


-.592442517 


-45.5528224 


45.5555933 


-1.53132531 


-90.6319756 


-.592415237 


-37.8144398 


37.3290272 


-1.53493114 


-93.7411959 


-.502375444 


-31.393595 


31.2974144 


-1.5367975 


-93.9143321 


-.592313242 


-24.0485212 


24.3453413 


-1.590 0 4995 


-91.1342735 


-.502235005 


-21.4322935 


21.4331279 


-1.59403911 


-9 1.3399271 


-.502114241 


-17.955175 


17.9421944 


-1.59375394 


-91.4818375 


-.501939134 


-14.9814731 


14.9999294 


-1.6344475 


-91.9292483 


-.531435399 


-12.3451825 


12.3753556 


-1.41134454 


-92.32339 


-.501317454 


-13.2532291 


10.2794714 


-1.41942733 


-92.7973*34 


-.530735075 


-3.59744911 


3.52219552 


-1.62959249 


-93.3433171 


-.533014484 


-7.05212434 


7.36932032 


-1.44153954 


-94.3554715 


-.493999591 


-5.84149977 


5.36295593 


-1.45599333 


-94.3314485 


-.497293211 


-4.33417139 


4.35943214 


-1.67338421 


-95.3734131 


-.494979437 


-3.99444734 


4.92519744 


-1.49407349 


-97.3435939 


-.491441952 


-3.29407913 


3.33054374 


-1.71395341 


-93.4891384 


-.434924784 


-2.78832265 


2.75174422 


-1.74363424 


-130.192243 


-.430213974 


-2.21734474 


2.26374904 


-1.73407334 


-102.219939 


-.478793143 


-1.30459822 


1.34580917 


-1.32599409 


-104.421937 


-.457774387 


-1.45455107 


1.52479424 


-1.37530954 


-137.447342 


-.448114424 


-1.14234732 


1.24233996 


-1.93275343 


-118.733943 


-.4147333 


-.913559365 


1.33414347 


-1.99831444 


-1 14.323672 


-.387925344 


-.793944749 


.832241914 


-2.37344959 


-113.381099 


-.359493245 


-.529223338 


.434372995 


-2.15631713 


-123.520724 


-.333454399 


-.334325031 


.494434454 


-2.24489756 


-123.422232 


-.242923931 


-.272339121 


.373992923 


-2.33771424 


-123.941293 


-.214412497 


-.184052021 


.235394043 


-2.43143359 


-139.314834 


-.172211374 


-. 122546394 


.211363185 


-2.5231131 


-144.56407 


-.132345421 


-.3731417391 


.154122443 


-2.49933337 


-149.535354 


-.3997494-31 


-.0434983292 


.113335129 


-2.43931891 


-154.115275 


-.9732317413 


-.32924S5321 


.0739939423 


-2.74134144 


-153.241914 


-.3529073594 


-.3172925105 


.3554424339 


-2.3254786 


-161.399516 


-.3377903073 


-.3190243947 


.3393112744 


-2.38166231 


-165.197174 


-.3244059837 


-5.78075S74E-93 


.3272393748 


-2. '9335 1322 


-147. 936334 


-.9134452453 


-3. 17315952E-03 


.313912329 


-2.97317314 


-179.250621 


-.0130013543 


-1.7121 1774E-33 


.3131134327 


-3.01965344 


-172.498385 


-9.93394734E-03 


-3.34429075E-34 


9.3771742E-03 


-3.34393135 


-17*. 487347 


-4.24175823E-93 


-4.22027R39E-04 


4.27402232E-93 


-3.37413774 


-176.135131 


-4.23274499E-83 


- 1.7 13892 19E-04 


4.3341503E-93 


-3.1820752 


-177.72533 


-2.99443421E-93 


-3.34033234E-05 


2.99448503E-93 


-3.12370133 


-179.241475 


-2.34773485E-33 


2.75143739E-35 


2.847947 12E-93 


3. 12323431 


179.237647 


-1.42707452E-33 


5.78349751E-95 


1.42821389E-03 


3.1014442 


177.7113 


-9.34502339E-04 


4.49614004E-95 


9. 3677693 IE-94 


3.07343161 


174.139947 


-4.73999033E-34 


4.4338 1543E-35 


4.3223082E-04 


3.34343413 


174.373361 


-4.48218744E-04 


6.19204743E-35 


4.72287539E-94 


3.91810415 


172.44444 


-3.22S19434E-94 


5.53937949E-95 


3.27434967E-94 


2.97255435 


179.315923 


-2.22557794E-04 


4.73502247E-05 


2. 276434 14E-94 


2.92931525 


147.364188 


-1.53422024E-04 


4.3924 1491E-05 


1.5379472E-34 


2.SS834356 


.00 . 0o *oo9 


-1.85743475E-94 


3.447579 13E-35 


1 . 1 130754SE-84 


2.32478945 


141.348571 


-7.29133322E-05 


2.91341137E-95 


7.35379394E-05 


2.7403443 


153.135929 


-5.32459559E-05 


2.44556444E-95 


5.5S936882E-05 


2.43377684 


154.955574 




Figure 8 
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£' TINDE 



Figure 8 (Cont.) 



PROBLEM IDENTIFICATION - EX 1 W' -PLANE T=.00 1 EXTENDED 

XX BODE PLOT (AMPLITUDE) XX 




ABSCISSA -> COMMON LOG OF FREQUENCY 
ORDINATE -> COMMON LOG OF AMPLITUDE 

MINIMUM FREQUENCY SHOWN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOWN ON ABSCISSA = 1OQ0 RADIANS/SEC 
AMPLITUDE LIMITS OF BODE PLOT ARE +-100 DECIBELS 



Figure 8 
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Figure 8 (Cont.) 



PROBLEM IDENTIFICATION - EX 1 W' -PLANE T=.Ci0 1 EXTENDED 

XX BODE PLOT (PHASE) XX 




ABSCISSA -> COMMON LOG OF FREQUENCY 
ORDINATE -> PHASE (DEGREES) 

TIC MARKS SHOW MULTIPLES OF 90 DEGREES 
MINIMUM FREQUENCY SHOWN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOWN ON ABSCISSA = 1000 RADIANS/SEC 
MAXIMUM PHASE ON ORDINATE SCALE = ISO DEGREES 
MINIMUM PHASE ON ORDINATE SCALE = -180 DEGREES 



Figure 8 
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FREQUENCY RESPONSE 



PROBLEM IDENTIFICATION - EX 1 W -PLANE T= . 0 1 

* x * x x * x x * x * x * * * x * x * * x x x x x * * x x x x x x x x x x x x x x xxxxxxxx x 



GAIN=-2 . 03 125 135E-05 

NUMERATOR COEFFICIENTS - IN ASCENDING POWERS OF S 



-2400400 . 84 
11302.0042 
1 



NUMERATOR ROOTS ARE 
REAL PART 



IMAGINARY PART 



-12002.0042 0 

200 0 



DENOMINATOR COEFFICIENTS - IN ASCENDING POWERS OF S 
0 

9.99167493 

1 

DENOMINATOR ROOTS ARE 

REAL PART IMAGINARY PART 

0 0 

-9.99167493 0 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 



Figure 9 w' -Plane Frequency Response Example, T=.01 
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Figure 


9 (cont . ) PR0BLEM identification 


- EX i w' 


-PLANE T= 


RADIAN FREQ. 


REAL PART 


IMAGINARY PART 


MAGNITUDE 


PHASE! RAD) 


PHASE! DEG) 


. 197226722 


-.524939547 


-46.6245341 


46.6274892 


-1.58285477 


-98.6450936 


.123284674 


-.524928388 


-40.5588649 


40.5534623 


-1.58374864 


-90.7416369 


.141747416 


-.524894366 


-35.2665654 


35.2704713 


-1.5856789 


-90.3527412 


.162975084 


-.524860366 


-38.6709756 


30.6754662 


-1.53790731 


-90.9384198 


.187381743 


-.524815428 


-26.6736542 


26.6788168 


-1.59846926 


-91.1272086 


.215443469 


-.524756033 


-23.1966296 


23.2825644 


-1.5934146 


-91.2959641 


.247787636 


-.524677539 


-20.1720793 


20.1789816 


-1.5968906 


-91.4899677 


.284803588 


-.52457381 


-17.5418867 


17.5488488 


-1.63869303 


-91.7129879 


.327454917 


-.524436749 


-15.2528909 


15.2611045 


-1.60516741 


-91.9693511 


.376493582 


-.524255672 


-13.2686853 


13.2710444 


-1.61031838 


-92.2640214 


.43287613 


-.524016491 


-11.527947 


11.5398508 


-1.61622129 


-92.6826918 


.497702359 


-.523790643 


-10.8200796 


18.3337559 


-1.62301398 


-92.9918344 


.572236769 


-.523233694 


-8.78767398 


8.72338307 


-1.63081873 


-93.439064 


.657933228 


-.52273353 


-7.56513572 


7.58317406 


-1.63978443 


-93.9527686 


.756463332 


-.522008018 


-6.57018639 


6.59089881 


-1.65088885 


-94.3427823 


.869749008 


-.52105292 


-5.78343855 


5.72718215 


-1.66190113 


-95.2199546 


1.00000001 


-.519793686 


-4.9479794 


4.97523789 


-1.67546416 


-95.9970594 


1.149757 


-.51813936 


-4.2891243 


4.32038737 


-1.69101639 


-96.3831655 


1.32194116 


-.515963635 


-3.71405452 


3.74972327 


-1.73883614 


-97.9091336 


1.51991189 


-.513126826 


-3.21161357 


3.25234698 


-1.72922959 


-99.077593 


1.74752841 


-.509417817 


-2.77209186 


2.81858943 


-1.75253532 


-108.412913 


2.99923302 


-.584596244 


-2.38784633 


2.43979662 


-1.77911839 


-101.93604 


2.31012972 


-.49836076 


-2.34916171 


2.1888924 


-1.309367 


-183.66913 


2.65608731 


-.493358539 


-1.75212349 


1.31944598 


-1.34367592 


-185.634337 


3.85385554 


-.488 148485 


-1.49052834 


1.56595572 


-1.88243386 


-107.355598 


3.51119177 


-.467296383 


-1.25981191 


1.34368586 


-1.92593991 


-119.351133 


4.3370173 


-.451325898 


-1.35619372 


1.14858184 


-1.97462499 


-113.137719 


4.64153888 


-.431817282 


-.376623211 


.97721191 


-2.82849743 


-116.224333 


5.33669929 


-.408476439 


-.718746596 


.326718149 


-2.93759426 


-119.610383 


6.13590734 


-.331235751 


-.580778561 


.694719472 


-2.15167917 


-123.29218 


7.35483239 


-.350349774 


-.461388989 


.579324936 


-2.22825434 


-127.211277 


8.1113084 


-.316458228 


-.359537358 


.47897969 


-2.29255366 


-131.353696 


9.32603358 


-.238578464 


-.274266592 


.392368887 


-2.36756996 


-135.651815 


10.7226724 


-.244087236 


-.20446534 


.318348855 


-2.444138 


-140.838842 


12.3234676 


-.208143963 


-.148767675 


.255843176 


-2.5210461 
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Figure 9 (Cont.) 



PROBLEM IDENTIFICATION - EX1 W' -PLANE T= . 0 1 

XX BODE PLOT (AMPLITUDE) XX 




ABSCISSA -> COMMON LOG OF FREQUENCY 
ORDINATE -> COMMON LOG OF AMPLITUDE 

MINIMUM FREQUENCY SHOWN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOWN ON ABSCISSA = 100 RADIANS/SEC 
AMPLITUDE LIMITS OF BODE PLOT ARE +-60 DECIBELS 



Figure 9 
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Figure 9 (Cont.) 



PROBLEM IDENTIFICATION - EX 1 W' -PLANE T= . 0 1 

XX BODE PLOT (PHASE) XX 




ABSCISSA -> COMMON LOG OF FREQUENCY 
ORDINATE -> PHASE (DEGREES) 

TIC MARKS SHOW MULTIPLES OF 90 DEGREES 
MINIMUM FREQUENCY SHOWN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOWN ON ABSCISSA = 100 RADIANS/SEC 
MAXIMUM PHASE ON ORDINATE SCALE = ISO DEGREES 
MINIMUM PHASE ON ORDINATE SCALE = -183 DEGREES 
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FREQUENCY RESPONSE 
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Figure 10 w' -Plane Frequency Response Example, T=.l 
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Figure 10 (Cont. )rroblEM IDENTIFICATION - EX 1 W' - PLANE T= . 
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Figure 10 (Cont.) 



PROBLEM IDENTIFICATION - EX 1 W' -PLANE T= . 1 

XX BODE PLOT (AMPLITUDE) XX 




ABSCISSA -> COMMON LOG OF FREQUENCY 
ORDINATE -> COMMON LOG OF AMPLITUDE 

MINIMUM FREQUENCY SHOWN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOWN ON ABSCISSA = 100 RADIANS/SEC 
AMPLITUDE LIMITS OF BODE PLOT ARE +-40 DECIBELS 



Figure 10 
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Figure 10 (Cont.) 



PROBLEM IDENTIFICATION - EX 1 W'-PLANE T= . 1 

XX BODE PLOT (PHASE) XX 




ABSCISSA -> COMMON LOG OF FREQUENCY 
ORDINATE -> PHASE (DEGREES) 

TIC MARKS SHOW MULTIPLES OF 90 DEGREES 
MINIMUM FREQUENCY SHOWN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOWN ON ABSCISSA = 100 RADI ANS/SEC 
MAXIMUM PHASE ON ORDINATE SCALE = 180 DEGREES 
MINIMUM PHASE ON ORDINATE SCALE = -180 DEGREES 



Figure 10 



109 



ROOT LOCUS 

PROBLEM IDENTIFICATION - EX 1 S-PLANE 

NUMERATOR COEFFICIENTS - IN ASCENDING POWERS OF S 
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Figure 11 s-Plane Root Locus Example 
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Figure 11 (Cont.) 
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Figure 11 (Cont.) 
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Figure 11 (Cont.) 
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Figure 11 (Cont.) 



PROBLEM IDENTIFICATION - EX 1 S-PLANE 
XX ROOT LOCUS PLOT XX 




ABSCISSA -> REAL (SIGMA) AXIS 
ORDINATE -> IMAGINARY (OMEGA) AXIS 
TIC MARKS SHOW INTERVALS OF 1 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -22 TO 1 1 
ORDINATE, -17 TO 14 
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Figure 12 w* -Plane Root Locus Example, T=.001 
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Figure 12 (Cont.) 
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Figure 12 (Cont.) 
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Figure 12 (Cont.) 
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Figure 12 (Cont . ) 



PROBLEM IDENTIFICATION - EX 1 W' -PLANE T=.00 1 
XX ROOT LOCUS PLOT XX 




ABSCISSA -> REAL (SIGMA) AXIS 
ORDINATE -> IMAGINARY (OMEGA) AXIS 
TIC MARKS SHOW INTERVALS OF 1 
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ORDINATE, -17 TO 16 



Figure 12 
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Figure 13 (Cont.) 
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Figure 13 (Cont.) 



14 


GAIN 


SS 


-7.0544 1 162E-66 


ROOTS 


ARE 










REAL PART 


I MAG. PART 










-2. 1954399 


0 










-7.71304879 


0 


15 


GAIN 


= 


-8.3179305 IE-06 


ROOTS 


ARE 










REAL PART 


I MAG. PART 










-2.82447443 


0 










-7.06911459 


0 


16 


GAIN 


= 


-9.77097722E-06 


ROOTS 


ARE 










REAL PART 


I MAG. PART 










-3.97303464 


0 










-5.90341974 


0 


17 


GAIN 


=: 


-1 . 144 198 IE-05 


ROOTS 


ARE 










REAL PART 


I MAG. PART 










-4.92837473 


-1 .73235163 










-4.92837473 


1 .73235163 


13 


GAIN 


= 


-1 . 33636352E-05 


ROOTS 


ARE 










REAL PART 


I MAG. PART 










-4.91704436 


-2.3109046 










-4.91704436 


2.3109046 


19 


GAIN 


= 


- 1 . 55735377E-05 


ROOTS 


ARE 










REAL PART 


I MAG. PART 










-4.90401439 


-3.65156925 










-4.9040 1439 


3.65156925 


20 


GAIN 


= 


-1 .811 49255E-0 5 


ROOTS 


ARE 










REAL PART 


I MAG. PART 










-4.33902985 


-4.42507143 










-4.88902985 


4.42507143 


21 


GAIN 


=: 


-2. 103752 14E-05 


ROOTS 


ARE 










REAL PART 


I MAG . PART 










-4.37179752 


-5. 1735032 










-4.87179752 


5. 1735032 


22 


GAIN 


= 


-2. 439850 68E-05 


ROOTS 


ARE 










REAL PART 


I MAG. PART 










-4.85193023 


-5.91326942 










-4.85198023 


5.91826942 


23 


GAIN 


= 


-2 . 32636399E-05 


ROOTS 


ARE 










REAL PART 


I MAG. PART 










-4.829190 13 


-6.67269346 










-4.829190 13 


6.67269846 


24 


GAIN 


= 


-3 . 2708543E-05 


ROOTS 


ARE 










REAL PART 


I MAG. PART 










-4.80293141 


-7.4463113 










-4.80298141 


7.44631 13 


25 


GAIN 


= 


-3.7820 18 16E-05 


ROOTS 


ARE 










REAL PART 


I MAG. PART 










-4.77234103 


-8.24663677 










-4.77234103 


3.24663677 


26 


GAIN 




-4 . 3698566E-05 


ROOTS 


ARE 










REAL PART 


I MAG. PART 










-4.73817921 


-9.080 10553 








_. -4.73817921 

Figure 13 


9.030 10553 



122 



Figure 13 (Cont.) 
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Figure 13 (Cont.) 



PROBLEM IDENTIFICATION - EX1 W y -PLANE T= . 0 1 
XX ROOT LOCUS PLOT XX 




ABSCISSA -> REAL (SIGMA) AXIS 
ORDINATE -> IMAGINARY (OMEGA) AXIS 
TIC MARKS SHOW INTERVALS OF 1 
THE PLOT FRAME LIMITS ARE; 
ABSCISSA, -22 TO 1 1 
ORDINATE , -17 TO 14 



Figure 13 
124 



ROOT LOCUS 

PROBLEM IDENTIFICATION - EX 1 W' -PLANE T= . 1 

X*X***X*X*XX*XX* XX* ********** *X**XX*XSSXXX**XXXX** 



NUMERATOR COEFFICIENTS - IN ASCENDING POWERS OF S OR W' 

-2439.71742 
10 1.985371 
1 



OPEN-LOOP 2EROS 

REAL PART IMAGINARY PART 

-121.935871 9 

20 0 



DENOMINATOR COEFFICIENTS - IN ASCENDING POWERS OF S OR W' 
0 

9.24234314 

1 



OPEN-LOOP POLES 



REAL PART 



IMAGINARY PART 



0 

-9.24234314 



0 

0 



MIN GAIN 
0 



MAX GAIN 
-. 113643529 




Figure 14 w' -Plane Root Locus Example, T=.l 



125 



ROOTS ARE 



Figure 14 (Cont.) 
1 GAIN = 0 



2 



3 



4 



5 



6 



7 



3 



6 



10 



1 1 



12 



13 



GAIN 



GAIN 



GAIN 



GAIN 



GAIN 



GAIN 



GAIN 



GAIN 



GAIN 



GA I N 



GAIN 



GAIN 



-1 .S870329SE-04 



-4.0 1422S4E-04 



-6 .•43344583E-04 



-9.3230 4545E-0 4 



- 1 . 25385852E-03 



-1 . 8343958E-03 



-2 . 08626324E-03 



-2.56291 102E-03 



-3. 13405597E-03 



-3 . 79087267E-03 



-4.54821 138E-03 



-5.414351 94E-03 



REAL PART 
-9.24234314 
0 

ROOTS 

REAL PART 
-.0496547373 
-9. 17536913 

ROOTS 

REAL PART 
-. 107695744 
-9.09740308 

ROOTS 

REAL PART 
-. 175741487 
-9 . 0084329 1 

ROOTS 

REAL PART 

-.255807277 

-8.89993997 

ROOTS 

REAL PART 

-.350442074 

-3.77500293 

ROOTS 

REAL PART 

-.482932327 

-8.62753258 

ROOTS 

REAL PART 

-.597624754 

-8.45263907 

ROOTS 

REAL PART 

-.780451059 

-3.24358792 

ROOTS 

REAL PART 
-.959374079 
-7 . 990S9 134 

ROOTS 

REAL PART 
- 1 .20372442 
-7.830 70 196 

ROOTS 

REAL PART 
-1 .52831341 
-7.29047243 

ROOTS 

REAL PART 



-1 .95977078 
-6.77784531 
Figure 14 



I MAG . 
0 
0 

ARE 

I MAG. 
0 
0 

ARE 

I MAG. 
0 
0 

ARE 

I MAG. 
0 
0 

ARE 

I MAG. 
0 
0 

ARE 

I MAG. 
0 
0 

ARE 

I MAG. 
0 
0 

ARE 

I MAG. 
0 
0 

ARE 

I MAG. 
0 
0 

ARE 

I MAG. 
0 
0 

ARE 

I MAG. 
0 
0 

ARE 

I MAG . 
0 
0 

ARE 

I MAG. 
0 



PART 



PART 



PART 



PART 



PART 



PART 



PART 



PART 



PART 



PART 



PART 



PART 



PART 



126 



0 



Figure 14 (Cont.) 

14 GAIN = -6.4 1 37380 3E -03 



15 GAIN = -7 . 56256453E-03 



16 GAIN = -8 . 3836575E-03 



17 GAIN = -.0104929144 



18 GAIN = -.0 121500599 



19 GAIN = -.0 141592772 



ROOTS 

REAL PART 

-2.61028387 

-6.03338216 

ROOTS 

REAL PART 

-4.26780979 

-4.26780979 

ROOTS 

REAL PART 

-4.20552827 

-4.20552827 

ROOTS 

REAL PART 
-4 . 13369394 
-4 . 13369394 

ROOTS 

REAL PART 

-4.05082203 

-4.05082208 

ROOTS 

REAL PART 



20 



21 



22 



23 



24 



25 



26 



GAIN 



GAIN 



GAIN 



GAIN 



GAIN 



GAIN 



GAIN 



-.0 16469377 



-.0 191270669 



-.022 1828352 



-.0256969633 



-.0297332224 



-.034385664 



-3.95515054 

-3.95515054 

ROOTS 

REAL PART 

-3.34464502 

-3.84464502 

ROOTS 

REAL PART 
-3.71692006 
-3.7 1692006 

ROOTS 

REAL PART 
-3.56917313 
-3.569 17313 

ROOTS 

REAL PART 

-3.39812932 

-3.39312932 

ROOTS 

REAL PART 
-3. 19939 14 
-3. 1998914 

ROOTS 

REAL PART 
-2.96936645 



-.039730222 



-2.96986645 

ROOTS 

REAL PART 
-2.70258524 
-2.70258524 
Figure 14 



127 



ARE 

I MAG. PART 
0 
0 

ARE 

I MAG. PART 
-.613935262 
• 6 1 39 352*2 

ARE 

I MAG . PART 

-2.04485029 

2.04435029 

ARE 

I MAG. PART 
-2.92566367 
2 . 92566367 

ARE 

I MAG. PART 
-3.6375661 1 
3 > o 3 / 5 £> 1 1 

ARE 

I MAG. PART 

-4.40426744 

4.40426744 

ARE 

I MAG. PART 
-5. 10621413 
5. 10621413 

ARE 

I MAG. PART 
-5.8102537 
5.9102587 

ARE 

I MAG. PART 

-6.52752509 

6.52752509 

ARE 

I MAG. PART 
-7.26633187 
7 . 26633 187 

ARE 

I MAG. PART 

-3.03350963 

8.03350963 

ARE 

I MAG. PART 

-3.33507658 

3.83507658 

ARE 

I MAG. PART 
-9.67661833 
9 . 67 66 1 888 



Figure 14 


(Cont . ) 






27 


GAIN 


sr 


-.0458764635 


ROOTS 


ARE 










REAL PART 


I MAG. PART 










-2.39151005 


-10. 5635 ! 88 










-2.39151005 


18 .5635 138 


28 


GAIN 


= 


-.0529446414 


ROOTS 


ARE 










REAL PART 


I MAG. PART 










-2.02878202 


-1 1.50 10951 










-2.02878202 


1 1 .50 10 951 


2? 


GAIN 


= 


-.0610730459 


ROOTS 


ARE 










REAL PART 


I MAG. PART 










-1 .60489341 


-12.4946347 










-1 .6048934 1 


12.4946347 


30 


GAIN 


= 


-.070420711 


ROOTS 


ARE 










REAL PART 


I MAG. PART 










-1 . 10825704 


-13.549675 










-1 . 10825704 


1 3 . 549675 


31 


GAIN 


=r 


-.081 170526 


ROOTS 


ARE 










REAL PART 


I MAG. PART 










-.52463290 1 


-14.6714325 










-.52463290 1 


14.6714325 


32 


GAIN 


= 


-.0935328132 


ROOTS 


ARE 










REAL PART 


I MAG. PART 










. 163647556 


-15.3654531 










. 163647556 


1 5 . 36545S 1 


33 


GAIN 


= 


-. 107749443 


ROOTS 


ARE 










REAL PART 


I MAG. PART 










.97874845 


-17. 1366797 










.97874845 


17. 1366797 



Figure 14 



128 



Figure 14 (Cont.) 
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C. ROOT LOCUS TEMPLATES IN w' -PLANE 



The s-plane consists of an abscissa representing the 
real portion of the complex variable s and an ordinate 
representing the imaginary part where 

s= 0" + jC 0 

The w' -plane approaches the s-plane as the sampling 
period T approaches zero. (See [Ref. 11] for a proof.) 

The w' -plane, like the s-plane, consists of an abscissa 
representing the real part and an ordinate representing 
the imaginary part where 

w'=u'+j V' 

The s and w' planes are related as shown below. 

u' = 2/T tanh ( <J T/2) 

V' = 2/T tan ( QJ ' T/2) 

<7 = 2/T tanh -1 (u'T/2) where -1 < u'T/2 < 1 

COj = 2/5 tan -1 ( ' T/2) where -1 < J^'T/2 < 1 

In the s-plane, (J represents damping. Therefore a line 
of constant damping in the s-plane is a vertical line 
parallel to the imaginary axis and passing through the proper 
value of O on the real axis. Since the real axis in the 
w'-plane represents u' which is not the damping but only 

related to C T as shown above, it is not possible to easily 

determine the damping of a particular root by simply ob- 
serving the position of the root in the plane. This is also 
true of the damped natural frequency, the natural frequency, 
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and the damping ratio. Because of this it is helpful to 
have templates of constant parameters in the w' -plane to 
make interpretation of the characteristics of a root 
by its location easier. The following sections give some 
insight into the nature of these templates in the w' -plane. 

1 . Constant Damping 

The constant damping templates are shown in Figure 15. 
There are three templates, one of each of three values of 
sampling period, .001 seconds, .01 seconds, and .1 seconds. 
Each template shows lines of constant damping for the values 
25, 50, 75, and 100. These templates are created by using 
the above relationship between (J and u' to find the constant 
value of u' that corresponds to the constant value of (J and 
then plotting the constant u' value. From the templates it 
is seen that for a period of .001 there is practically no 
difference from the s-plane. For a period of .01 seconds, 
some distortion is noticed as the values of damping approach 
the reciprocal of the value of the period. This distortion 
is in the form of u' becoming smaller than the value of 
damping. And for a period of .1 seconds gross distortion 
is seen. Figure 15 shows an extra template for a period 
of .1 seconds but for values of damping of 1, 2, 3, and 4. 

In this case only little distortion is present because the 
values of damping are small compared to the reciprocal of 
the period. 
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2 . Constant Damped Natural Frequency 



Figure 16 shows the templates for constant damped 
natural frequency for values of damped natural frequency of 
25, 50, 75, and 100. The results here are similar to the 
constant damping case. As the period increases, the 
distortion between the value of OJj and the value of P* in- 
creases. In this case the value of P l becomes greater than 
the value of damped natural frequency as it distorts. 

In the case where the period is .1 seconds the distortion 
is so great that the template is useless and negative 
values of P 1 are produced by the tangent function in the 
relationship . 

3 . Constant Natural Frequency 

Templates for constant natural frequency in the 
w' -plane are shown in Figure 17. In the s-plane constant 
natural frequency CJ n plots as a circle with its center at 
the origin where 

OJh = sqrt ( CT 2 + CJj 2 ) 

To plot constant natural frequency in the w' -plane, substitute 



u* = 2/T tanhC <7 T/2) 
P' = 2/T tan ( T/2) 



into the above equation and solve for P 1 to get 



P' = 2/T tan 



T/2 (CJ 2 - [2/T tanh 1 (u'T/2)] 2 }^ 
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For a constant value of natural frequency and period this 
equation is used to plot u' vs y* . In the w' -plane for a 
period of .001 seconds and a value of natural frequency of 25 
it plots very close to a perfect circle. As the values of 
natural frequency increase they form concentric ellipses 
which are elongated along the imaginary axis and contracted 
along the real axis. This effect becomes more dramatic 
as the period increases. An extra template for a period of 
.05 seconds is added to better show the trend since the 
change from .01 seconds to .1 seconds is so great. 

4 . Constant Damping Ratio 



formed by the negative real axis and a line joining the 
origin and the root in a root locus plot. In the s-plane a 
plot of constant damping ratio is a line radially out from 



real axis. To plot constant damping ratios in the w' -plane, 
substitute the proper relations into the above equation to 
get 



The damping ratio, f , can be defined as 



£ = ( CT / GJn | = cos Q 




the origin forming the correct angle $ with the negative 



£ 



tanh 1 (u'T/2) 



{[tanh 1 (u'T/2)] 2 + [tan 1 (^'T/2)] 2 p 
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Solving for u' gives 



u' = 2/T tanh 



tan 1 ( T/2 ) £ 

sqrt ( 1- £ 2 ) 



This equation is used to plot constant damping ratio for 
a given sampling period. Figure 18 shows templates for 
constant damping ratio for values of damping ratio of .1, 

.5, .707, and .9. From these templates it can be seen that 

there is a region near the origin that can be interpreted 
the same as the s-plane for a close approximation and that 
the size of this region depends on the value of the sampling 
period. Beyond this region a template makes the interpreta- 
tion easier. 
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CONSTANT DAMPING W' -PLANE T=.801 

XX ROOT LOCUS PLOT XX 




ABSCISSA -> REAL <U> AXIS 
ORDINATE -> IMAGINARY (NU> AXIS 
TIC MARKS SHOW INTERVALS OF 5 
THE PLOT FRAME LIMITS ARE : 
ABSCISSA, -163 TO 37 
ORDINATE, -109 TO 100 



Figure 15 Constant Damping Templates 



Figure 15 (Cont.) 



CONSTANT DAMPING W' -PLANE T= . 0 1 

XX ROOT LOCUS PLOT XX 




ABSCISSA -> REAL ( U) AXIS 
ORDINATE -> IMAGINARY CNU) AXIS 
TIC MARKS SHOW INTERVALS OF 5 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -169 TO 40 
ORDINATE, -100 TO 100 



Figure 15 
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Figure 15 (Cont.) 



CONSTANT DAMPING W'-PLANE T= . 1 

XX ROOT LOCUS PLOT XX 




ABSCISSA -> REAL < U) AXIS 
ORDINATE -> IMAGINARY (NU) AXIS 
TIC MARKS SHOW INTERVALS OF 5 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -118 TO 82 
ORDINATE, -100 TO 100 



Figure 15 
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Figure 15 



CONSTANT DAMPING W y -PLANE T= . 1 

XX ROOT LOCUS PLOT XX 




ABSCISSA -> REAL ( U) AXIS 
ORDINATE -> IMAGINARY <NU) AXIS 
TIC MARKS SHOW INTERVALS OF 1 
THE PLOT FRAME LIMITS ARE; 
ABSCISSA, -7 TO 3 
ORDINATE, -5 TO 5 



Figure 15 
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CONSTANT DAMPED NATURAL FREQ W' -PLANE T=.001 
XX ROOT LOCUS PLOT XX 




ABSCISSA -> REAL (U) AXIS 
ORDINATE -> IMAGINARY <NU) AX I 
TIC MARKS SHOW INTERVALS OF 5 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -100 TO 180 
ORDINATE, -37 TO 163 



S 



Figure 16 Constant Damped Natural Frequency Templates 
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Figure 16 (Cont.) 



CONSTANT DAMPED NATURAL FREQ W' -PLANE T= . 0 1 
XX ROOT LOCUS PLOT XX 



J ITT IT 1 1 1 1 1 1 1 II 1 1 1 1 1 


T'mTTT'T 111)11 ITT ITT 


XJ.l.l 11 11 t l II 111 11 II 1 


i i i i i i i i it i i i i i i ( m r 



ABSCISSA -> REAL <U> AXIS 
ORDINATE -> IMAGINARY (NU) AXIS 
TIC MARKS SHOW INTERVALS OF 5 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -100 TO 100 
ORDINATE, -32 TO 163 



Figure 16 
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Figure 16 (Cont.) 



CONSTANT DAMPED NATURAL FREQ W' -PLANE T= . 1 
XX ROOT LOCUS PLOT XX 



_l 1 III' M M 1 1 111 ill 1 1 
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ABSCISSA -> REAL <U> AXIS 
ORDINATE -> IMAGINARY (NU) AXIS 
TIC MARKS SHOW INTERVALS OF 5 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -100 TO 100 
ORDINATE, -105 TO 95 



Figure 16 
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CONSTANT NATURAL FREQUENCY W / -PLANE T=.001 
XX ROOT LOCUS PLOT XX 




ABSCISSA -> REAL <U> AXIS 
ORDINATE -> IMAGINARY (NU) AXIS 
TIC MARKS SHOW INTERVALS CF 5 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -153 TO 52 
ORDINATE, -105 TO 100 



Figure 17 Constant Natural Frequency Templates 
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Figure 17 (Cont.) 



CONSTANT NATURAL FREQUENCY W' -PLANE T= . 0 1 
XX ROOT LOCUS PLOT XX 




ABSCISSA -> REAL (U) AXIS 
ORDINATE -> IMAGINARY <NU> AXIS 
TIC MARKS SHOW INTERNALS OF 5 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, - 15G TO 55 
ORDINATE, -105 TO 108 



Figure 17 
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Figure 17 (Cont.) 



CONSTANT NATURAL FREQUENCY W' -PLANE T=.05 
XX ROOT LOCUS PLOT XX 




ABSCISSA -> REAL ( U) AXIS 
ORDINATE -> IMAGINARY <NU) AXIS 
TIC MARKS SHOW INTERVALS OF 5 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -128 TO 77 
ORDINATE, -105 TO 100 



Figure 17 
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Figure 17 (Cont.) 



CONSTANT NATURAL FREQUENCY W' -PLANE T= . 1 
XX ROOT LOCUS PLOT XX 




ABSCISSA -> REAL < U) AXIS 
ORDINATE -> IMAGINARY CNU> AXIS 
TIC MARKS SHOW INTERVALS OF 5 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -120 TO 85 
ORDINATE, -105 TO 100 



Figure 17 
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CONSTANT DAMPING RATIO W' -PLANE T=.801 
XX ROOT LOCUS PLOT XX 




ABSCISSA -> REAL < U) AXIS 
ORDINATE -> IMAGINARY <NU) AXIS 
TIC MARKS SHOW INTERVALS OF 100 
THE PLOT FRAME LIMITS ARE: 

ABSCISSA, -3000 TO 1800 
ORDINATE, 0 TO 4000 



Figure 18 Constant Damping Ratio Templates 
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Figure 18 (Cont.) 



CONSTANT DAMPING RATIO W' -PLANE T= . 0 1 
** ROOT LOCUS PLOT ** 




ABSCISSA -> REAL <U> AXIS 
ORDINATE -> IMAGINARY <NU) AXIS 
TIC MARKS SHOW INTERVALS OF 10 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -360 TO 100 
ORDINATE, 0 TO 400 



Figure 18 
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Figure 18 (Cont.) 



CONSTANT DAMPING RATIO W' -PLANE T= . 1 
XX ROOT LOCUS PLOT XX 
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ABSCISSA -> REAL (U) AXIS 
ORDINATE -> IMAGINARY <NU> AXIS 
TIC MARKS SHOW INTERVALS OF 1 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -30 TO 10 
ORDINATE, 9 TO 40 



Figure 18 
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D. COMPENSATION EXAMPLE 



To further demonstrate the use of the RTLOC and FRESP 
programs for analysis and design in the w' -plane, an example 
from [Ref. 6] is duplicated here. 

An angular position servo is described below. 




where 



K = 1.9 (V/rad) 

P 

K t = 0.385 (Nm/A) 



K t = 1(A/V) 

J = 3.22 x 10 3 (Nm. Sec 2 ) 



Neglecting friction (Md) gives 




where 



ZOH 



1-e 



-sT 



s 
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G 1 (s) = 



KjK t 1-e 



-sT 



Converting to the z domain gives 



G 1 Cz) = K 2 T" 



( z+1 ) 
(z-1 )‘ 



where 



K 



2 



KjK t 



2 J 



59.8 



Converting this to the w domain gives 
K,T 2 1— w 

G i ( »> - 4- — 

2 w 

where 

w = w'T/2 

The open loop transfer function is 
G(w) = K K 1 G 1 (w) 

p 1 1 

where 

K x = 3.83 
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from the calculation of D C gain. This gives 



G(w ) = 217.58 T 



1 — w 



w 



For 



T = 0.01 

1-w 

G(w ) = 0.0218 — (5) 
w 



Therefore for the uncompensated system equation 5 above is 
entered into the FRESP program producing the output shown 
in Figure 19. From these Bode plots it is seen that the 
uncompensated system is unstable. 

From analog design techniques the final compensation 
network [Ref. 7] is 

(w/0.1 + 1) 

K-D(w) = 3.83 

(w/2.15 + 1)“ 

and the open loop transfer function becomes 

w+0.1 w-1 

G(w) = K p — «- 

(w+2.15) w 



For 



K = -1.017 



151 



this transfer function was entered into the FRESP program. 

The resulting Bode plots in Figure 20 show that the compen- 
sated system is stable with a gain margin of 10.5 dB and 
phase margin of 41 degrees. 

Figure 21 shows the output from the RTLOC program for the 
compensated system. This figure includes the normal root 
locus in addition to an expanded portion of the root locus 
to obtain more detail. These plots are of w and not w' in 
this case. Values of w' can be found by the simple relation 

w' = 2w/T 
where 

T = .01 seconds 
for this problem. 

On the expanded root locus plot in Figure 21, a root 
associated with a certain gain may be identified and 
templates similar to the ones discussed in section V.C. may 
be used to easily obtain the characteristics associated with 
the root. 
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FREQUENCY RESPONSE 

PROBLEM IDENTIFICATION - KATZ EXAMPLE 

XXXXXXXX X X X X X X X X X X X X X XXX X X X X X X X X X X X X X X X XXX XXX X X X X 7 

GAIN=- . 022 

NUMERATOR COEFFICIENTS - IN ASCENDING POWERS OF S 



-1 

1 



NUMERATOR ROOTS ARE 
' REAL PART 



1 



IMAGINARY PART 
0 



DENOMINATOR COEFFICIENTS - IN ASCENDING POWERS OF 

9 

0 

1 



DENOMINATOR ROOTS ARE 
REAL PART 



0 

9 



IMAGINARY PART 

0 

0 



¥ ¥ ¥ X ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ ¥ 

A A /\ A A A A ✓ A /•. /\ /s \ A ✓ *. /\ ✓ % A ✓s /\ /\ A /\ , \ .*% /*. A . \ A* ,-s / \ \ 



Figure 19 Uncompensated System Frequency Response 
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Figure 19 (Cont.) 

PROBLEM IDENTIFICATION - KATZ EXAMPLE 



RADIAN FREQ. 


REAL PART 


IMAGINARY PART 


MAGNITUDE 


PHASE (RAD) 


PHASE(DEG) 


.184761575 


-2.88455806 


.21088066 


2.31552385 


3.83721133 


174.819434 


.1149757 


-1.66421932 


.19134473 


1.6751832 


3.82711966 


173.441242 


.126135688 


-1.33166411 


. 174346238 


1.39262867 


3.91687843 


172.303163 


.138488637 


-1.14783182 


. 158857798 


1.15382953 


3.30397933 


172.115492 


.151991189 


-.952327477 


.144745309 


.963264672 


2.99875683 


171.357759 


.166318854 


-.798639801 


.131336535 


.80 1563527 


2.97628452 


178.529743 


.183873329 


-.656482388 


.128178898 


.667311733 


2.968524 


169.625591 


.288923381 


-.544956794 


.109494513 


.555347962 


2.94339951 


163.639273 


.228513875 


-.452432784 


.8997673267 


.463382133 


2.92455386 


167.564607 


.242812328 


-.375617579 


.0989842724 


.336461862 


2.90414539 


166.295333 


.26568878 


-.311844312 


.33232S5873 


.322656861 


2.33197322 


1 05.125243 


.291585393 


-.253393627 


.8754783239 


.269674375 


2.35794731 


163.743377 


.319926716 


-.214942188 


.9637657432 


.225674261 


2.33195624 


162.259193 


.351119176 


-.17844373 


.3626567336 


.139129163 


2.39392119 


163.652908 


.335352362 


-.148151312 


.3570905323 


.153778716 


2.77377659 


153.925743 


.422924291 


-.122997322 


.3523137667 


.133545559 


2.74148151 


157.375376 


.464153387 


-.192114953 


.3473975628 


.112573323 


2.78792683 


155.181263 


.599413386 


-.8347776276 


.3431363939 


.8951433593 


2.67844244 


153.335136 


.559881823 


-.8703838757 


.8393582893 


.8386370536 


2.63139424 


158.791329 


.613598733 


-.8534339381 


.835854513 


.0635578613 


2.5912399 


143.467163 


.673415872 


-.3435123753 


.8326693815 


.3534374545 


2.54393266 


146.043126 


.739872211 


-.3482762413 


.3297670507 


.8500324612 


2.38512214 


143.532977 


.31113884 


-.8334338431 


.827122628 


.9438558773 


2.46818142 


148.953479 


.398215895 


-.3277608588 


.3247131234 


.9371672377 


2.41421 


123.324393 


.977909969 


-.323847546 


.8225176322 


.3322216693 


2.3673227 


135.666296 


1.87226724 


-.9191344776 


.9285172734 


.0239553662 


2.32133516 


133.392755 


1.17681197 


-.3153857794 


.0136945753 


.3245325322 


2.27514727 


128.356383 


1.29154968 


-.3131336531 


.3178338803 


.32154277 


2.22964625 


127.749266 


1.41747413 


-.8109494515 


.3155285649 


.8189941631 


2. 185199S3 


125.202259 


1.55567617 


-9.898 42782E-83 


.3141417637 


.3163114621 


2. 14289359 


122.733252 


1.70735267 


-7.54783221E-83 


.3123854456 


.3149329292 


2 . .986ooo5 


128.35764 


1.37381745 


-6.26567372E-93 


.3117407331 


.31330383 


2.06131571 


113.337544 


2.85651234 


-5.28137654E-83 


.3106977233 


.311395411 


2.32339131 


115.931352 


2.25791976 


-4.31363929E-03 


9.74736784E-93 


.0136612495 


1.93736531 


113.396362 


2.4770764 


-3.58545172E-03 


3.33143732E-93 


9.57735998E-83 


1.95448931 


111.934857 


2.71853329 


-2.9767858 IE-03 


3.39243536E-93 


3.62254583E-03 


1.92327327 


113.195431 


2.93364729 


-2.47131279E-83 


7.3735257E-93 


7.77664891E-83 


1.39419926 


183.529146 


3.27454922 


-2.35172729E-03 


6.71848199E-93 


7.32473364E-93 


1 .36713663 


186.931955 


3.59381373 


-1.78338804E-93 


6.12163857E-03 


6.35420053E-83 


1.34213783 


135.549532 


3.94428613 


-1.41417683E-83 


5.57788 179E-93 


5.75428246E-83 


1.31918063 


134.226329 


4.32376137 


-1.17407373E-93 


5.93223524E-93 


5.21613578E-83 


1.79732662 


183.007915 


4.75081026 


-9.74736632E-04 


4.63878902E-03 


4.7322635SE-03 


1.77325322 


181.236727 


5.2140884 


-8. 892435 16E-84 


4.2 1940249E-03 


4. 2963845 IE-93 


1.76823637 


188.357816 


5.72236773 


-6.71843134E-94 


3.S445624E-83 


3.99232462E-83 


1.74338223 


99.9125428 


6.23829153 


-5.57738 164E-94 


3.50332287E-83 


3.5471513E-03 


1.72369989 


99.847197 


6.39261226 


-4. 6397839 IE-04 


3.19132324E-83 


3.22524372E-03 


1.71437396 


93.2558756 


7.56463344 


-3.34456231E-04 


2.98827946E-83 


2.93357133E-03 


1.70222342 


97.5385339 


3.30217537 


-3. 19132317E-94 


2.64999773E-83 


2. 66986 132E-83 


1.69866924 


96.3632463 


9.11162777 


-2.64998767E-84 


2.41449723E-93 


2.42899595E-83 


1.63818376 


96.2631754 


18.3888002 


-2. 19999992E-84 


2. 19999995E-03 


2.21097261E-83 


1.67046583 


95.71863 






Figure 19 
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Figure 19 (Cont.) 



PROBLEM IDENTIFICATION - KATZ EXAMPLE 

XX BODE PLOT (AMPLITUDE) XX 




ABSCISSA -> COMMON LOG OF FREQUENCY 
ORDINATE -> COMMON LOG OF AMPLITUDE 

MINIMUM FREQUENCY SHOWN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOWN ON ABSCISSA = 10 RADIANS/SEC 
AMPLITUDE LIMITS OF BODE PLOT ARE + -49 DECIBELS 



Figure 19 
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Figure 19 (Cont.) 



PROBLEM IDENTIFICATION - KATZ EXAMPLE 

XX BODE PLOT (PHASE) 




ABSCISSA -> COMMON LOG OF FREQUENCY 
ORDINATE -> PHASE (DEGREES) 

TIC MARKS SHOW MULTIPLES OF 90 DEGREES 
MINIMUM FREQUENCY SHOWN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOWN ON ABSCISSA = 10 RADIANS/SEC 
MAXIMUM PHASE ON ORDINATE SCALE = 180 DEGREES 
MINIMUM PHASE ON ORDINATE SCALE = 90 DEGREES 



Figure 19 
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FREQ U EN C Y R E S P ON S E 
PROBLEM IDENTIFICATION 

Z % z z z zzzzz z z z z z z z z z z z z z 



KATZ EXAMPLE 

zzzzzzzzzzzzzz 



XX XX 



GAIN=- 1 .017 



NUMERATOR COEFFICIENTS - IN ASCENDING POWERS OF S 
1 

-.? 

1 



NUMERATOR ROOTS ARE 
REAL PART 



IMAGINARY PART 



1 0 
1 9 



D EN OM I NAT 0 R COEFFICIENTS 



IN ASCENDING POWERS OF S 



0 

0 

4 . 6225 

4.3 

l 



DENOMINATOR ROOTS ARE 
REAL PART 



IMAGINARY PART 



- 2.15 0 

- 2.15 0 

0 0 

0 0 



XX X XXX X X X XXX X XX v XXXX X X X X X X xxxxxxxxxx X X X X X X X X X X X X X x 
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Figure 20 Compensated System Frequency Response 
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KATZ EXAMPLE 



Figure 20 (Coat.) 

PROBLEM IDENTIFICATION - 



RADIAN FREQ. 


REAL PART 


IMAGINARY PART 


MAGNITUDE 


PHASE (RAD) 


PHASE(CEG) 


. 184741575 


-2.39228771 


-1.66891537 


2.91226675 


-2.53478893 


-145.227717 


.1149757 


-2.85136854 


-1.58789532 


2.54547622 


-2.58797881 


-143.696143 


.126185688 


-1.76321733 


-1.366342S3 


2.23461116 


-2.43373669 


-142.333962 


. 133488637 


-1.5329249 


-1.23744316 


1.97885693 


-2.46245977 


-141.333337 


.151991189 


-1.33732181 


-1.11923564 


1.74391164 


-2.44471727 


-149.372332 


.166818854 


-1.17461536 


-1.81835445 


1.54969324 


-2.43398665 


-139.2S5325 


. 183873829 


-1.83916282 


-.911228484 


1.33289273 


-2.42169891 


-133.753176 


.283923381 


-.926259956 


-.319533617 


1.23676643 


-2.41725149 


-133.499353 


.220513875 


-.831990732 


-.735817878 


1.11816156 


-2.41388832 


-133.541263 


.242812828 


-.7536395 


-.656961126 


.999378661 


-2.42426292 


-133.988983 


.26568878 


-.686823804 


-.584713337 


.992018563 


-2.43632263 


-139.591957 


.291585308 


-.638984245 


-.517732353 


.316122474 


-2.45443313 


-149.62371 


.319926716 


-.533485553 


-.455338223 


.749839985 


-2.47382133 


-142.026851 


.351119176 


-.542781741 


-.397275793 


.672572193 


-2.53968973 


-143.794631 


.335352362 


-.587412466 


-.342966767 


.612443866 


-2.54721596 


-145.944776 


.422924291 


-.476353936 


-.292885357 


.553777222 


-2.59155943 


-143.434955 


.464158387 


-.44852322 


-.244322736 


.518755483 


-2.64231151 


-151. 421999 


.509413386 


-.42304587 


-.199427817 


.46769463 


-2.79137383 


-154.760472 


.559881823 


-.399158742 


-.157212169 


.428995316 
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Figure 20 
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Figure 20 (Cont.) 



PROBLEM IDENTIFICATION - KATZ EXAMPLE 

XX BODE PLOT (AMPLITUDE) XX 




ABSCISSA -> COMMON LOG OF FREQUENCY 
ORDINATE -> COMMON LOG OF AMPLITUDE 

MINIMUM FREQUENCY SHOWN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOWN ON ABSCISSA = 10 RADIANS/SEC 
AMPLITUDE LIMITS OF BODE PLOT ARE +-60 DECIBELS 



Figure 20 



Figure 20 (Cont.) 



PROBLEM IDENTIFICATION - KATZ EXAMPLE 

XX BODE PLOT (PHASE) XX 




ABSCISSA -> COMMON LOG OF FREQUENCY 
ORDINATE -> PHASE (DEGREES) 

TIC MARKS SHON MULTIPLES OF 90 DEGREES 
MINIMUM FREQUENCY SHOWN ON ABSCISSA = .1 RADIANS/SEC 
MAXIMUM FREQUENCY SHOWN ON ABSCISSA = 10 RADIANS/SEC 
MAXIMUM PHASE ON ORDINATE SCALE = 180 DEGREES 
MINIMUM PHASE ON ORDINATE SCALE = -180 DEGREES 



Figure 20 
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ROOT LOCUS 

PROBLEM IDENTIFICATION - 

************ * * * * * * * * * * * * ; 


KATZ EXAMPLE 

(*** * * * * * * * * * * * * * * * * * * * * * 


NUMERATOR COEFFICIENTS - 
1 

_ Q 

1 


IN ASCENDING POWERS OF 3 


OPEN-LOOP ZEROS 
REAL PART 


IMAGINARY PART 


1 

1 


0 

0 


D EN OM 1 NAT OR C 0 E F F I C I ENT 3 

0 

0 

4.6225 
4 . 3 
1 


- IN ASCENDING POWERS OF 


OPEN-LOOP POLES 
REAL PART 


IMAGINARY PART 


-2. 15 
-2.15 
0 
0 


0 

0 

0 

0 


MIN GAIN 
0 


MAX GAIN 
-30 


* * * * * * * * * * * 6 * * ******** * * r 

OPTIQN TAKEN 


C * * * * * * * * * * * * * * * * * * * * * * * * 


SIGMA MIN = -1 
OMEGA MIN = -1 


SIGMA MAX = 1 
OMEGA MAX = 1 



Figure 21 Compensated System Root Locus 
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Figure 21 (Cont.) 



PROBLEM IDENTIFICATION - KATZ EXAMPLE 
XX ROOT LOCUS PLOT XX 



f 1 I I 1 I I 

/ 


1 


\ 

■ i i i i i i 


I ‘' 1 mhh>+ +++ ■ 
1 



ABSCISSA -> REAL (SIGMA) AXIS 
ORDINATE -> IMAGINARY (OMEGA) AXIS 
TIC MARKS SHOW INTERVALS OF 1 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, -8 TO 2 
ORDINATE, -5 TO 5 



Figure 21 
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Figure 21 (Cont.) 



PROBLEM IDENTIFICATION - KAT2 EXAMPLE 

ROOT LOCUS PLOT XX 



— i 

J X 


— 1 


♦ 

♦ 

\1 


's 



ABSCISSA -> REAL (SIGMA) AXIS 
ORDINATE -> IMAGINARY (OMEGA) AXIS 
TIC MARKS SHOW INTERVALS OF 1 
THE PLOT FRAME LIMITS ARE: 
ABSCISSA, - 1 TO 2 
ORDINATE, -1 TO 2 



Figure 21 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



This thesis demonstrates that common classical control 
programs may be adapted to run on an inexpensive microcomputer 
system. With the exception of the memory management problem 
encountered, this was done with relative ease. This allows 
more people access to these programs by introducing a new 
group of computers on which they can be run. Also these 
programs require no knowledge of any computer language or 
input card formats making them easier to use than the 
Fortran versions. 

It is also demonstrated that the transfer function pro- 
grams are a useful tool in the study of sampled data systems 
as well as classical systems. 

Only a sampling of five programs was converted in this 
thesis. Other programs can and should be similarly converted 
to run on microcomputer systems. There are programs in 
existence for the same microcomputer system used in this 
thesis that generate the aircraft stability derivatives. 

An effort should be made to modify these programs and the 
programs of this thesis to make them compatible and 
complementary . 
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APPENDIX A 



DESCRIPTION OF MICROCOMPUTER SYSTEM 

The microcomputer system used in developing this thesis 
consisted of the following components: 

Apple II plus computer (48K) 

Disk II 5 1/4" floppy disk drive and controller card 

USI 9" green screen monitor 

NEC PC-8023A-C dot matrix printer 

Grappler printer interface card 

Add Ram 16K expansion card 

The programming language used was Applesoft basic. All 
graphs were generated using the High Resolution graphics 
commands. In the High Resolution graphics mode a 
matrix of dots 280 dots wide and 192 dots high can be dis- 
played. The High Resolution page one is the only page used 
in this thesis and it resides in memory in the 8,192-byte area 
from $2000 to $3FFF. High Resolution page two resides in 
the area from $4000 to $6000. 

A memory management problem was encountered during the 
programming of the programs requiring the use of the graphics 
capabilities of the microcomputer system. To understand the 
nature of the problem encountered a brief description of the 
normal use of memory in the Apple II is necessary. An Apple- 
soft program is normally loaded at memory location $800 and 
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loads up. LOMEM is set to the end of the program. Simple 
variables are stored from LOMEM upward as they are defined 
in the program. This gives only 6K-bytes of RAM before there 
is a conflict with the first High-Resolution page and only 14K 
before there is a conflict with the second High-Resolution 
page. [Ref. 8] Since the programs and variables in this 
thesis exceed 14K, simply using High-Resolution page two in- 
stead of one is not the answer. The first part of the solu- 
tion was to get the Applesoft programs to load above the 
space used by High-Resolution page one and in effect pro- 
tecting that space from interference. This fix creates a 
new problem. The disk operating system loads at the top of 
the 48K of memory and sets HIMEM to $9CF8. String variables 
start at HIMEM and build down. This allows 22K of space for 
the program and all variables. Some programs in this thesis 
require more space than this. Even with a 16K memory 
expansion card installed the disk operating system will 
ignore it so the added memory is useless. This final problem 
was solved by using a utility program that relocated the disk 
operating system into the higher memory provided by the 
memory expansion card and resets HIMEM to $BFOO. This gives 
32. 5K of useable program space which is sufficient for all 
programs in this thesis. [Ref. 9] 
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